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ABSTRACT 

Steam-condensation  heat-transfer  measurements  were  made 
using  a  5-tube  in-line  test  condenser  with  an  additional 
perforated  tube  to  simulate  up  to  30  active  tubes.   Results 
were  obtained  for  smooth  tubes,  wire-wrapped  tubes  and 
dropwise-coated  tubes.   The  average  outside  heat-transfer 
coefficient  for  30  smooth  tubes  was  0.64  times  the  Nusselt 
coefficient  for  the  first  tube.   A  total  of  eight  wire- 
diameter  and  wire-pitch  combinations  were  tested:   1.6-mm- 
diameter  wire  wrapped  at  16  mm,  7.6  mm  and  4  mm  wire  pitches, 
1 . 0-mm-diameter  wire  wrapped  at  8  mm,  4  mm  and  2  mm  wire 
pitches,  and  0 . 5-mm-diameter  wire  wrapped  at  4  mm  and  2  mm 
wire  pitches.   The  best  bundle  performance  was  obtained  when 
the  tubes  were  wrapped  with  1 .0-mm-diameter  wire  at  a  wire 
pitch  of  4  mm.   This  combination  resulted  in  an  average  out- 
side heat-transfer  coefficient  for  30  tubes  that  was  1.15 
times  the  value  computed  for  the  first  tube  using  the  Nusselt 
theory.   The  average  outside  heat-transfer  coefficient  for 
the  30  dropwise-coated  tubes  was  1.1  times  the  value  of  the 
heat- transfer  coefficient  for  the  first  tube  in  the  tube 
bundle.   Utilizing  either  wire-wrapped  tubes  or  dropwise- 
coated  tubes,  it  is  possible  to  significantly  reduce  the 
condenser  surface  area  and  overall  size. 
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I.   INTRODUCTION 

A.   HISTORICAL  BACKGROUND 

Considerable  interest  has  been  generated  in  reducing  the 
size  and  the  weight  of  propulsion  systems  for  naval  applica- 
tions.  Advances  in  condenser  design  could  do  much  to  reduce 
the  size  and  the  weight  of  the  propulsion  plant.   Measures 
to  raise  the  condensing-side  heat-transfer  coefficient  (of 
condenser  tubes)  is  one  way  to  achieve  this  reduction  in 
condenser  size.   This  reduction,  however,  comes  at  a  price. 
This  is  usually  due  to  an  increase  in  the  pumping  power  or 
due  to  an  increase  in  the  initial  cost  of  the  tubes.   For 
naval  applications,  where  the  size  of  a  vessel  may  depend 
upon  the  size  of  the  condenser  (in  a  submarine,  for  example), 
this  reduction  in  the  size  is  justified,  even  at  the  greater 
cost. 

Search  [Ref.  1],  at  the  Naval  Postgraduate  School, 
investigated  the  present  condenser  design  process,  and 
examined  the  potential  benefits  that  might  occur  if  heat- 
transfer  enhancement  was  established  in  the  condenser.   He 
concluded  that  reductions  of  as  much  as  forty  percent  in  the 
size  and  weight  of  condensers  are  possible.   This  is  depen- 
dent, of  course,  on  the  heat-transfer-enhancement  technique 
utilized.   Much  further  research  work  at  the  Naval  Post- 
graduate School  has  been  directed  toward  these  heat-transfer- 
enhancement  techniques. 
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Beck  [Ref .  2] ,  Pence  [Ref .  3] ,  Reilly  [Ref .  4] ,  Fenner 
[Ref .  5] ,  and  Ciftci  [Ref.  6]  conducted  research  employing 
a  single- tube  test  condenser.   Their  research  concluded  that 
the  overall  heat-transfer  coefficient  of  enhanced  tubes  may 
be  as  much  as  twice  that  for  smooth  tubes  of  similr  geometry. 
In  a  separate  report,  Marto,  Reilly  and  Fenner  [Ref.  7]  re- 
ported that  most  of  the  increase  in  the  overall  heat-transfer 
coefficient  was  on  the  cooling-water  side  and  was  due  to  an 
increase  of  the  turbulence  and  the  swirl,  as  well  as  to  an 
increase  in  the  inside  surface  area.   Only  a  small  increase 
occurred  on  the  steam  side. 

Present-day  steam  condensers  utilizing  smooth  tubes  are 
limited  in  their  thermal  efficiency,  due  primarily  to  the 
large  thermal  resistances  occurring  on  the  tube  side  of  the 
condenser.   It  is  possible,  however,  by  utilizing  enhanced 
tubes,  to  increase  the  inside  heat- transfer  coefficient  by 
100  percent  or  more.   The  corresponding  increase  in  the  out- 
side heat-transfer  coefficient,  however,  is  less  than  50 
percent.   In  studying  ways  to  further  increase  the  outside 
heat-transfer  coefficient,  Webb  [Ref.  8]  reported  that 
conduction  across  the  condensate  film  is  the  primary  thermal 
resistance  in  film  condensation.   This  thermal  resistance  is 
usually  larger  than  the  thermal  resistance  of  the  tube  wall, 
that  attributed  to  fouling  or  that  due  to  noncondensable 
gases.   It  is  possible  to  reduce  this  thermal  resistance  by 
utilizing  a  geometry  that  reduces  the  film  thickness.   This 
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reduction  of  the  thermal  resistance  would  mean  a  correspond- 
ing increase  in  the  outside  heat-transfer  coefficient. 

For  large  tube  bundles,  condensate  inundation  is  present 
and  must  be  considered  when  attempting  to  increase  the 
overall  heat-transfer  coefficient.   Thomas  [Ref.  9]  wrapped 
wire  around  smooth  tubes  in  a  helical  manner  and  tested  the 
condensation  of  ammonia  in  a  large  tube  bundle.   Increases 
of  as  much  as  200  percent  in  the  oustide  heat-transfer 
coefficient  over  that  predicted  by  Nusselt  [Ref.  10]  were 
measured.   This  increase  was  attributed  to  the  effect  of 
surface  tension  drawing  the  condensate  to  the  wire  and 
acting  as  a  condensate  run-off  channel .   In  a  paper  by 
Cunningham  [Ref.  11],  "roped"  tubes  were  considered  and 
again  the  increase  in  the  outside  heat-transfer  coefficient 
was  attributed  to  condensate  drainage,  while  an  increase  in 
the  inside  heat-transfer  coefficient  was  attributed  to  the 
increase  in  the  inside  convective  coefficient,  due  to  the 
increased  turbulence.   Kanakis  [Ref.  12]  tested  both  smooth 
and  roped  tubes,  with  and  without  a  wire  wrap  in  an  in-line 
tube  bundle  simulating  up  to  30  tubes.   Adding  the  wire  wrap 
on  the  smooth  tube  increased  the  average,  outside  heat- 
transfer  coefficient  for  the  30  tube  bundle  by  50  percent, 
while  adding  the  wire  warp  to  the  roped  tube  increased  the 
average  outside  heat-transfer  coefficient  for  the  30  tube 
bundle  by  more  than  3  5  percent. 

For  condenser  tubes,  the  increase  in  the  outside  heat- 
transfer  coefficient  may  also  be  accomplished  by  promoting 
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dropwise  condensation.   Tanasawa  [Ref.  13]  reviewed  dropwise 
condensation  and  discussed  the  methods  for  promoting  these 
dropwise  conditions.   He  concluded  that  the  use  of  organic 
polymers  (such  as  Teflon)  was  the  most  promising  of  all  the 
dropwise-promoting  techniques.   Investigations  of  organic 
coatings  by  Brown  [Ref.  14]  have  shown  that  enhancements  of 
up  to  180  percent  are  possible.   Of  primary  concern,  however, 
is  that  these  coatings  have  very  low  conductivities  and  must 
therefore  be  ultra-thin.   Another  concern  is  that  these 
coatings  must  be  strongly  adherent  and  sufficiently  tough  to 
withstand  the  conditions  of  their  assembly  and  their  use. 
Holden  [Ref.  15]  investigated  numerous  dropwise-promoting 
coatings.   His  tests  were  based   not  only  on  the  ability  of 
the  coating  to  promote  dropwise  condensation,  but  also  on 
the  ability  of  the  coating  to  sustain  this  dropwise  perfor- 
mance over  an  extended  period  of  time.   In  addition,  the 
coatings  that  were  able  to  sustain  this  dropwise  performance 
were  evaluated  for  their  heat-transfer  performance.   His 
results  indicated  that  the  outside  heat-transfer  coefficient 
for  a  single  tube  could  be  increased  by  a  factor  of  from 
five  to  eight  through  the  use  of  polymer  coatings. 

Since  the  use  of  a  wire  wrap  can  increase  the  heat- 
transfer  performance  of  a  tube  bundle,  an  important  question 
arises:   Is  there  some  optimal  wire  diameter  and  wire  pitch 
combination?   In  addition,  during  dropwise  conditions  what 
is  the  effect  of  condensate  inundation  upon  heat  transfer  in 
a  tube  bundle? 
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B.   OBJECTIVES 

The  objectives  of  this  thesis  were  therefore  to: 

1.  Conduct  steam-condensation  tests  to  determine  the 
steam-side  heat-transfer  coefficient  for  16-mm-o.d. 
smooth  titanium  tubes, 

2.  Confirm  the  heat-transfer-performance  measurements 
of  Kanakis  [Ref.  12]  on  16-mm-o.d.  smooth  titanium 
tubes  with  a  1.6-mm-o.d.  titanium  wire  wrapped  at  a 
pitch  of  7.6  mm, 

3.  Conduct  heat-transfer-performance  measurements  to 
determine  the  effect  of  wire  pitch  and  wire 
diameter  on  the  steam-side  heat-transfer 
coefficient  with  inundation,  and 

4.  Conduct  heat-transfer-performance  measurements  to 
determine  the  effect  of  inundation  on  a  dropwise 
coated  tube  and  compare  this  performance  with  the 
performance  of  a  smooth  tube. 
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II.   THEORETICAL  AND  EMPIRICAL  BACKGROUND 

The  basis  for  the  analysis  of  film  condensation  on  a 
horizontal  tube  was  set  forth  by  Nusselt  in  1916.   His  analy- 
sis was,  however,  for  laminar  film  condensation  on  a  single 
horizontal  tube.   Nusselt' s  analysis  yielded  the  well-known 
relationship  for  the  heat-transfer  coefficient: 


h„   =   0.725 
Nu 


V3   2  , 

kf  Pf  hfq  g 

Do  ^f(Tsat-V 


1/4  (2.1) 


This  relationship  is  subject  to  the  following  restric- 
tions, as  stated  by  Nobbs  [Ref.  16]: 

1.  The  wall  temperature  is  constant, 

2.  The  flow  is  laminar  in  the  condensate  film, 

3.  The  film  thickness  is  small   compared  to  normal 
tube  diameters, 

4.  All  fluid  properties  are  constant  within  the 
condensate  film, 

5.  Heat  transfer  in  the  film  is  by  conduction, 

6.  All  forces  due  to  hydrostatic  pressure,  surface 
tension,  inertia,  and  vapor/liquid  interfacial  shear 
are  negligible  when  compared  to  viscous  and 
gravitational  forces,  and 

7.  The  vapor/liquid  interface  and  the  surrounding  steam 
are  at  the  saturation  temperature. 

Jakob  [Ref.  17]  extended  the  Nusselt  analysis  to  film 

condensation  on  a  vertical  in-line  column  of  horizontal 

tubes  by  assuming  that  all  the  condensate  from  a  tube  drains 
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as  a  laminar  sheet  onto  the  tube  below  it.   This  is  depicted 
in  Figure  2.1a.   In  this  idealized  situation,  the  average 
coefficient  for  a  vertical  column  of  N  tubes  was  predicted 
to  be : 


h>7   =   0.725 

N 


kf  Pf  hfq  g 

Do  N  *r  (Tsat'V 


1/4 

(2.2) 


Combining  equations  (2.1)  and  (2.2)  yields  the  Nusselt 
idealized  theory  for  the  average  coefficient  compared  to  the 
coefficient  of  the  top  tube: 


hN/hx   =   N  1/4  (2.3) 


In  terms  of  the  local  heat-transfer  coefficient  for  the  Nth 
tube,  this  result  becomes: 


hN/h1   =   N3/4  -  (N  -1)3/4  (2.4) 


Realizing  that  condensate  does  not  flow  in  a  laminar  sheet 
but  by  discrete  droplets,  Kern  [Ref.  18]  proposed  a  less- 
conservative  relationship  to  account  for  the  ripples  and 
turbulence  introduced  into  the  condensate  film.   This  is 
depicted  in  Figure  2.1b.   The  Kern  relationship  is: 


h  /h,   =   N  1/6  (2.5) 

N   1 
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a.   Nusselt's  Idealized 
Laminar  Flow  Model 


b.   Kern's  More  Realistic 
Flow  Model 


c.   Eissenberg's  Side-Drainage  Flow  Model 


Figure  2.1   Representations  of  Condensate  Fl 


ow 
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or  in  terms  of  the  local  heat-transfer  coefficient  for  the 
Nth  tube: 


hN/h1   =   N5/6  -  (N  -1)5/6  (2.6) 


Eissenberg  [Ref.  19]  did  extensive  experimentation  to 
investigate  the  effects  of  steam  velocity,  condensate  inun- 
dation and  noncondensable  gases  on  the  condensation  heat- 
transfer  coefficient.   He  theorized  that  condensate  does  not 
always  drain  only  onto  tubes  aligned  vertically,  but  can  be 
diverted  sideways,  especially  in  staggered  tube  bundles. 
The  condensate  draining  sideways  strikes  the  lower  tubes  on 
their  sides  rather  than  on  their  tops.   This  is  depicted  in 
Figure  2.1c.   Since  more  heat  is  transferred  from  the  top  of 
a  tube  than  from  its  bottom,  the  net  effect  of  inundation  is 
less  severe.   He  predicted  the  following  relationship: 


hN/hl   =   °*60  +  °'42  N  1/4  (2'7) 


Extensive  experimental  research  into  the  effect  of 
condensate  inundation  has  been  conducted.   However,  the  data 
exhibit  a  substantial  amount  of  scatter  as  shown  in  Figure 
2.2  by  the  cross-hatched  area.   Berman  [Ref.  20]  conducted  a 
compilation  of  film  condensation  data  on  bundles  of  horizontal 
tubes  and  identified  the  following  variables  as  important 
in  causing  the  scatter: 
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1.  Bundle  geometry  (in-line  or  staggered), 

2.  Tube  spacing, 

3.  Type  of  condensing  fluid, 

4.  Operating  pressure, 

5.  Heat  flux,  and 

6.  Local  vapor  velocity. 

Marto  and  Wanniarachchi  [Ref.  21]  noted  that  other  factors 
can  have  an  effect  on  the  data  such  as  noncondensable  gases, 
the  direction  of  the  vapor  flow,  partial  dropwise  condi- 
tions, an  insufficient  abount  of  steam  reaching  the  lower 
tubes  in  the  bundle  and  difficulty  in  the  measurement  of 
the  local  condensing  coefficient. 

Very  small  amounts  of  noncondensable  gases  can  result  in 
significant  reductions  in  the  condensation  heat-transfer 
rate.   Summaries  of  the  phenomenon  have  been  made  by  Chisholm 
[Ref.  22],  and  Webb  and  Wanniarchchi  [Ref.  23].   The  consen- 
sus among  them  is  that  noncondensable  gases  impose  an  added 
thermal  resistance,  since  the  vapor  molecules  must  diffuse 
through  a  gas  layer  prior  to  reaching  the  condensing  surface. 
These  noncondensable  gases  can  also  lead  to  regions  where 
the  tubes  are  inoperative  in  a  condensing  role.   Nonconden- 
sable gases  have  an  adverse  effect  on  the  condenser  perfor- 
mance, and  must  be  taken  into  account  when  designing  a  tube 
bundle . 

The  motion  of  vapor  within  a  condenser  affects  the  film 
condensation  process  because  of  its  effect  on  the  surface 
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shear  between  the  vapor  and  the  film,  and  the  resulting  effect 
on  vapor  separation.   Although  the  results  may  be  different 
depending  upon  the  orientation  of  the  steam  flow,  the  general 
effect  of  an  increase  in  vapor  velocity  is  a  corresponding 
increase  in  the  condensing  heat-transfer  coefficient.   There 
exist  a  number  of  both  theoretical  and  empirical  equations 
representing  vapor-shear  effects  on  the  condensate  heat- 
transfer  coefficient  [Refs.  24,  25  and  26]. 

It  is  clear  that  the  measurements  that  were  made  during 
this  thesis  had  both  inundation  and  vapor-shear  effects. 
Marto  [Ref.  27]  stated  that  these  two  effects  are  difficult 
to  separate  from  one  another.   During  this  study,  however, 
an  attempt  was  made  to  separate  these  two  effects.   For  this 
purpose,  data  were  taken  on  the  top  tube  of  the  bundle  with 
vapor  velocity  as  a  variable,  and  the  data  were  correlated 
using  an  expression  similar  to  a  correlation  suggested  by 
Fujii  [Ref.  28]. 


1/2 

Nu/Re~A     =   BF  (2.8) 

2(p 


This  correlation  is  based  on  numerous  data  for  both 
staggered  and  in-line  tube  bundles.   To  represent  the  test- 
condenser  tube  bundle,  the  constants  B  and  n  in  equation 
(2.8)  were  computed  using  a  least-squares  technique. 
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III.   EXPERIMENTAL  APPARATUS 

The  test  facility  was  designed  and  built  originally  by 
Morrison  [Ref .  29] ,  and  modified  by  Noftz  [Ref .  30]  and 
Kanakis  [Ref.  12]  to  simulate  an  active  tube  column  of  up  to 
30  in-line  tubes.   A  detailed  description  of  all  the  com- 
ponents is  contained  in  Ref.  30  and  further  modifications 
are  described  in  Ref.  12.   The  descriptions  in  this  report 
are  brief,  with  particular  focus  given  to  the  modifications 
undertaken  by  the  author. 

A.  STEAM  SUPPLY 

The  steam  supply  system  is  shown  in  Figure  3.1.   House 
steam  flows  through  a  supply  valve  (MS-3) ,  into  a  steam 
separator,  and  is  throttled  down  by  another  valve  (MS-4) 
to  operating  conditions.   The  steam  then  passes  through  an 
orifice  and  into  the  test  condenser  diffuser  before  entering 
the  test  condenser  where  it  flows  through  a  simulated  in-line 
tube  bundle . 

B.  TEST  CONDENSER 

The  dimensions  of  the  test  condenser  shown  in  Figure  3.2 
were  unchanged  from  Noftz' s  original  design.   The  condenser 
consisted  of  five  active  tubes,  twelve  dummy  tubes  and  one 
perforated  tube.   The  tubes  were  positioned  in  the  test  con- 
denser with  a  pitch-to-diameter  ratio  of  1.5.   The  active 
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Fiaure  3.1   Schematic  Diagram  of  Test  Apparatus 


26 


^-79wmfr| 


TEST 

SECTION 


EXHAUST 
PLENUM 


PERFORATED 
TUBE 


ACTIVE 
TUBES 


NOTE1   ALL  COMPONENTS  DRAWN  TO  SCALE 


Figure  3.2   Sketch  of  Test  Condenser 
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condensing  length  of  each  tube  was  305  mm.   A  double-walled 
glass  viewing  window  was  provided  on  the  front  of  the  test 
condenser  to  allow  for  visual  observation  of  the  condensation 
process.   Although  the  viewing  window  was  designed  to  have 
heated  air  flowing  between  the  glass  walls,  this  feature  was 
not  utilized  since  no  fogging  of  the  glass  was  observed.   A 
modification  of  the  test  condenser  was  made  to  minimize  the 
steam  flow  along  the  test  condenser  side  walls.   The  modifica- 
tion consisted  of  the  addition  of  twelve  baffle  plates,  placed 
as  illustrated  in  Figure  3.2. 

C.   TEST-CONDENSER  TUBES 

Ten  different  types  of  active  tubes  were  tested  in  this 
set  of  experiments.   The  tubes  were  manufactured  by  Wolverine 
Division  of  Universal  Oil  Products,  and  were  made  of  titanium. 
All  were  of  16  mm  o.d.  and  had  a  minimum  wall  thickness  of 
0.89  mm.   The  ten  types  of  tubes  are  as  listed  below: 

1 .  Smooth  tubes , 

2.  Smooth  tubes  wrapped  with  1.6-mm-o.d.  titanium  wire 
at  three  different  pitches  (16  mm,  7.6  mm  and  4  mm), 

3.  Smooth  tubes  wrapped  with  1.0-mm-o.d.  titanium  wire 
at  three  different  pitches  (8  mm,  6  mm  and  4  mm) , 

4.  Smooth  tubes  wrapped  with  0.5-mm-o.d.  titanium  wire 
at  two  different  pitches  (4  mm  and  2  mm) ,  and 

5.  Smooth  tubes  with  a  dropwise  coating. 

The  actual  wrapping  of  the  wire  around  the  smooth  tubes 
was  performed  by  the  author  of  this  thesis.   A  detailed 
description  of  the  procedure  is  included  in  Chapter  III. 
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The  dropwise-coated  tubes  were  sent  to  General  Magna- 
plate  Corporation  where  a  dropwise  coating,  commercially 
referred  to  as  "Nedox,"  was  applied.   First  an  electro- 
deposited  nickel-cobalt  substrate  was  applied  to  the  tube 
surface,  and  then  this  substrate  was  infused  with  a  micro- 
thin  Teflon  layer.   A  heat-treating  cycle  was  employed  to 
ensure  thorough  infusion  of  the  substrate  with  the  Teflon. 
This  coating  is  highly  non-wetting,  has  a  low  coefficient 
of  friction  and  is  heat-resistant.   The  specifics  of  the 
"Nedox"  process  are  proprietory,  but  the  manufacturer  claims 
control  of  the  surface  thickness  to  0.0025  mm.   For  this  thesis 
the  requested  surface  thickness  was  0.0075  mm. 

D.   PERFORATED  TUBE 

Each  set  of  five  active  tubes  had  its  own  corresponding 
perforated  tube.   This  perforated  tube  was  located  above  the 
uppermost  active  tube;  see  Figure  3.2.   Figure  3.3  shows  a 
typical  pattern  for  the  perforations  of  a  tube,  and  is  an 
illustration  of  the  supply-water  tubing  at  the  test  condenser. 
This  particular  arrangement  of  the  tubing  was  used  in  an 
attempt  to  provide  an  even  flow  from  the  entire  length  of  the 
perforated  tube.   A  heated  tube  was  used  to  provide  the 
water  supply  for  the  perforated  tube.   A  separate  rotameter 
was  provided  in  order  to  control  the  rate  of  flow  of  condensate 
through  the  perforated  tube. 
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E.  CONDENSATE  SYSTEM 

The  condensate  system  remained  unchanged  from  Kanakis ' 
modifications.   Figure  3.1  shows  the  system  as  designed  by 
Noftz  and  as  modified  by  Kanakis.   As  steam  condenses  in  the 
test  condenser,  it  runs  into  and  collects  in  the  test  con- 
denser hotwell.   The  hotwell  is  a  calibrated  cylinder  equipped 
with  a  sight  glass,  providing  a  visual  measure  of  the 
condensate  being  produced.   Measurement  of  the  condensate 
collection  rate  provided  the  flow  rate  to  be  supplied  to 
the  perforated  tube. 

F.  COOLING  WATER  SYSTEM 

Although  the  cooling  water  system  remained  unchanged 
from  Noftz  original  design,  the  five  rotameters  were  recali- 
brated after  Kanakis'  data  were  taken.   Each  active  tube  had 
its  own  rotameter  and  regulating  valve,  and  it  was  possible 
to  control  the  coolant  velocity  through  the  active  tubes  up 
to  values  of  5  m/s.   All  data  for  this  thesis,  however,  were 
taken  at  a  constant  coolant  velocity  of  1.56  m/s.   The  inlet 
temperature  of  the  cooling  water  was  maintained  at  a  nearly 
constant  temperature  by  the  use  of  a  large  supply  tank  and 
a  separate  cooling  tower  to  dissipate  excess  heat  to  the 
atmosphere.   In  order  to  provide  a  bulk  temperature  at  the 
outlet  of  each  active  tube,  a  mixing  chamber  was  installed 
just  upstream  of  the  temperature  monitoring  point. 
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G .   INSTRUMENTAT ION 

1 .  Flow  Rates 

Rotameters  were  used  to  measure  the  flow  rate  of  the 
cooling  water  to  each  active  tube  and  to  the  perforated  tube. 

2 .  Temperatures 

Stainless-steel-sheathed  copper-constantan  thermo- 
couples were  used  as  the  primary  temperature  monitoring 
devices.   Table  1  lists  the  locations  monitored  by  each  of 
the  thermocouples.   Calibration,  as  performed  by  Kanakis, 
was  utilized  for  all  data  taken. 

The  perforated-tube  water  supply  was  not  maintainable 
at  a  constant  temperature,  nor  was  it  quickly  adjustable  to 
a  desired  temperature.   The  system,  as  designed,  contained 
a  Gulton  Industries,  West  20  temperature  controller.   Although 
the  manufacturer's  accuracy  was  listed  as  1.25  degrees  F, 
the  system  response  was  so  poor  that  the  temperature  of  the 
perforated-tube  supply  water  varied  by  as  much  as  10  degrees 
C  from  the  beginning  of  the  data  run  to  the  completion  of 
the  data  run. 

3 .  Pressure 

A  Bourdon-tube  pressure  gauge  was  used  to  measure  the 
house  steam  supply  pressure,  while  a  compound  pressure  gauge 
was  used  to  measure  the  steam  pressure  downstream  of  the 
orifice.   A  pressure  transducer  was  used  to  measure  the 
pressure  in  the  test  condenser,  and  a  U-tube  mercury  manom- 
eter was  fitted  to  measure  the  pressure  drop  across  the 
orifice. 
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TABLE  1 
Thermocouple  Monitoring  Locations 
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4 .   Data  Collection  and  Display 

A  Hewlett-Packard  model  3054A  Automatic,  Data- 
Acquisition  System,  with  a  HP  2671G  printer  was  used  to 
record  and  display  the  thermocouple  and  transducer  readings 
The  pressure  transducer  was  assigned  to  channel  19  of  the 
data-acquisition  system,  with  the  thermocouples  assigned  to 
the  channels  as  indicated  in  Table  1. 
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IV.   EXPERIMENTAL  PROCEDURES 

A.  CALIBRATION 

1.  Rotameters 

Calibration  was  performed  on  the  five  rotameters 
that  supplied  the  cooling  water  to  the  active  tubes.   This 
calibration  was  performed  using  a  weighing  tank  and  a  stop- 
watch, and  a  least-squares-fit  calibration  line  was  generated 
for  each  of  the  rotameters. 

2 .  Steam  Orifice 

To  calibrate  the  steam  orifice  the  following  steps 
were  taken: 

1.  Set  a  steam  flow  rate, 

2.  Adjust  the  cooling  water  flow  rate  to  condense  all 
the  incoming  steam, 

3.  Measure  the  pressure  drop  across  the  orifice, 

4.  Measure  the  condensate  collection  rate, 

5.  Repeat  steps  1  through  4  for  different  steam  flow 
rates,  and 

6.  Develop  a  least-squares-fit  straight  line  for  the 
steam  mass  flow  rate  vs.  pressure  drop. 

B.  PREPARATION  OF  CONDENSER  TUBES 

Prior  to  the  winding  of  the  titanium  wire,  the  outside 
surface  of  the  tubes  was  buffed  using  steel  wool  unitl  all 
evidence  of  surface  oxidation  was  removed.   If  any  oxidation 
was  evident  on  the  inner  surface  of  the  tube,  a  test-tube 
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brush  was  used  to  apply  a  50-percent  sulfuric-acid  solution 
to  the  inner  surface  of  the  tube.   The  tube  was  then 
thoroughly  rinsed  using  tap  water. 

The  procedure  followed  in  the  wrapping  of  the  wire  onto 
the  smooth  tubes  was  the  same  for  every  wire  diameter  and 
pitch  combination  tested.   The  wire  was  first  welded  to  the 
tube  and  then,  under  a  constant  tension  of  22  Newtons ,  the 
wire  was  guided  (using  a  metal  plate  to  control  the  spacing 
between  the  wire  wraps)  to  the  proper  pitch  as  the  tube  was 
manually  turned.   Upon  completion  of  the  wrapping,  the  free 
end  of  the  wire  was  then  welded  to  the  tube  to  yield  a 
helically-wrapped  surface  305  mm  in  length.   After  the  wire 
wrapping  of  the  tubes  was  completed  and  just  prior  to  the 
installation  of  the  tubes  in  the  test  condenser,  the  tubes 
were  cleaned  inside  and  out  using  a  brush  and  a  biodegradable 
detergent.   A  thorough  rinsing  was  agin  undertaken,  and  the 
tubes  were  checked  to  ensure  that  they  displayed  proper  wetting 
characteristics.   If  any  non-wetting  areas  were  evident,  a 
50-percent  solution  of  sodium  hydroxide  mixed  with  an  equal 
amount  of  ethyl  alcohol  and  heated  to  about  80  degrees  C,  was 
brushed  onto  the  tube  surface.   After  rinsing  with  tap  water, 
the  tube  was  ready  for  installation  in  the  test  condenser. 
Upon  completion  of  the  installation  of  the  tubes  in  the  test 
condenser,  steam  was  introduced  into  the  test  condenser. 
Again  the  tubes  were  visually  examined.   If  there  were  any 
areas  that  displayed  dropwise  condensation,  the  tubes  were 
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washed  from  above  using  the  perforated- tube  water  supply. 
This  was  sufficient  in  all  cases  to  restore  filmwise 
condensation. 

C.   SYSTEM  OPERATION 

The  system  was  operated  in  accordance  with  the  operating 
instructions  outlined  in  Appendix  A  of  Ref.  17.   The  system 
was  deemed  to  be  operating  at  a  steady-state  condition  when 
the  cooling  water  inlet  temperature  was  steady.   The  cooling 
water  inlet  temperature  was  monitored  using  the  HP-3054A 
Automatic  Data  Acquisition  System,  and  was  considered  steady 
when  two  consecutive  readings  (at  five-minute  intervals) 
varied  by  less  than  0.1  degrees  C. 

One  hour  was  usually  required  from  system  start-up  to 
the  steady-state  condition.   Upon  reaching  the  steady-state 
condition,  a  data  run  was  begun.   Any  changes  in  the  cooling 
water  flow  or  changes  in  the  perforated-tube  water  flow  were 
accompanied  by  a  five-minute  wait  for  system  restabilization, 
prior  to  taking  any  further  data.   The  duration  of  each  data 
set  was  approximately  one  minute,  and  a  completed  data  run 
consisted  of  30  separate  data  sets.   During  each  data  set, 
the  condensation  collection  rate  was  measured.   For  each 
inundation  condition,  five  consecutive  data  sets  were  taken 
and  the  average  values  were  computed.   After  each  set  of 
five  consecutive  data  sets,  the  average  condensate  collection 
rate  was  used  to  determine  the  rate  of  flow  of  the  water  into 
the  perforated  tube  which  was  to  be  used  for  the  next  set  of 


37 


data  (the  average  amount  of  condensate  collected  during  a 
data  set  was  introduced  into  the  perforated  tube  for  the 
next  inundation  condition).   For  example,  to  simulate  tubes 
6  through  10,  the  perforated-tube  water  flow  rate  was  set 
equal  to  the  condensate  collection  rate  for  tubes  1  through 
5.   A  data  run  was  considered  completed  upon  the  simulation 
of  inundation  conditions  through  the  30th  tube. 

Thermocouple  readings  and  pressure-transducer  readings 
were  taken  automatically  by  the  data-acquisition  system, 
while  the  settings  of  the  rotameters  and  the  test  condenser 
initial  and  final  hotwell  levels  were  entered  into  the 
computer  using  the  keyboard.   Initial  and  final  hotwell 
levels  were  taken  for  a  time  interval  of  one  minute. 

The  following  conditions  were  the  operating  conditions 
for  all  data  taken  for  this  thesis: 

1.  The  coolant  velocity  was  1.55  ±  0.05  m/s, 

2.  The  coolant  inlet  temperature  was  24  degrees  C.  (±3 
degrees  C.  depending  on  the  day  of  the  run,  but 
constant  during  a  given  run) ,  and 

3.  The  saturation  temperature  was  100.5  ±  0.5  degrees  C. 

D.   DATA- REDUCTION  PROGRAM 

A  computer  program  was  utilized  to  process  and  plot  all 
raw  data.   The  program  was  written  in  BASIC  language  and 
was  run  on  an  HP-9  826  computer.   A  listing  of  this  program 
is  included  as  Appendix  C. 
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E.   HEAT-TRANSFER-COEFFICIENT  CALCULATION 

1.  Inside  Heat-Transfer  Coefficient 

The  inside  heat-transfer  coefficient  was  computed 
using  the  Sieder-Tate  equation  described  in  Holman  [Ref. 
31]  : 

Nu   =   h.  D./kc   =   Ci  Re0'8  Pr1/3  (yc/yw)°*14  (4.1) 

where  the  coefficient  C.  was  computed  using  a  modified 
Wilson  plot  [Ref.  12].   Using  this  technique,  Kanakis  [Ref. 
12]  obtained  a  Sieder-Tate  coefficient  of  0.029  ±  0.001. 
The  method  Kanakis  used,  however,  did  not  include  the  variation 
of  the  condensate  film  properties  as  a  function  of  heat  flux. 
A  slightly  different,  modified  Wilson  plot,  as  described  by 
Nobbs  [Ref.  16],  that  does  account  for  this  variation  in 
the  film  properties,  was  used  in  this  thesis.   Appendix  D 
gives  a  description  of  this  method.   Reprocessing  Kanakis1 
data  with  this  method,  the  Sieder-Tate  coefficient  was  deter- 
mined to  be  0.028  ±  0.001. 

2 .  Outside  Heat-Transfer  Coefficient 

The  heat-transfer  rate  to  the  cooling  water  can  be 
computed  using  an  energy  balance,  and  can  be  expressed  in 
terms  of  an  overall  heat-transfer  coefficient  by: 


Q   =   mc    (T    -T.)   =   U   A   LMTD  (4.2) 

C   pC    CO    CI         o   o 


where 
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mTD      "      ln((T  -=%    J77T   „-T   ))  (4-2a: 

sat   ci     sat   co 


Solving  for  the  overall  heat-transfer  coefficient  gives 


U    =   ra   c   /A   ln((T    -T  .)/(T    -T   ))  (4.3) 

o      c   pc/  o      sat   ci"   sat    co  '  v    ' 


The  overall  thermal  resistance  can  be  written  as  the  sum  of 
the  internal,  wall  and  external  resistances: 


+  ^  +  ^r~  (4.4) 


U   A       h.  A.    A     h   A 
O   O       11     o     o   o 


where 


R    =   D   ln(D  /D. )/2  k  (4.4a) 

w      o     o   l     m 


or  from  this  result,  the  outside  heat-transfer  coefficient 
may  be  written: 


ho   -   1/U   -  d'/D.  h.  -  R  <4-4b» 

O      Oil      w 


F.   DATA-REDUCTION  TECHNIQUE  FOR  THE  OUTSIDE  HEAT-TRANSFER 
COEFFICIENT 

A  listing  of  the  computer  program  used  in  the  calculation 

of  the  outside  heat-transfer  coefficient  can  be  found  in 

Appendix  C.   The  steps  required  in  the  technique  are  outlined 

below: 
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1.   Calculate  the  average  bulk  cooling  water  temperature: 


T,   =   (T    +  T  .  )/2 
b       co    ci  ' 


(4.5) 


2.   Calculate  the  cooling  water  velocity 


V    =   m  /A.  p 
c       c'  1  Kc 


(4.6) 


3.   Calculate  the  cooling  water  Reynolds  number 


Re   =   p   V  D./y 
Kc   c  v     c 


(4.7) 


4.  Calculate  the  heat  transferred  to  the  cooling 
water  using  the  left-hand  side  of  equation  (4.2) 

5.  Calculate  the  heat  flux: 


q   =   Q/(tt  Dq  L) 


(4.8) 


6.   Since  the  film  temperature  was  not  a  known  quantity, 
an  iterative  scheme  was  employed  to  calcualte  the 
film  temperature,  as  indicated  below: 

a)  Assume  a  film  temperature  (say  T   =  T   .), 

b)  Calculate  the  Nusselt  coefficient, 


Nu 


=   0.651 


Pkf  pf  hfg  j 

vf  °o  * 


1/3 


(4.9) 


c)   Evaluate  Tr   ,  and 
f  ,c 


f,c      sat    2  h 


Nu 


(4.10) 
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d)   Repeat  steps  a)  through  c)  until  T_    is 
approximately  equal  to  T  .  'c 

7.   Calculate  the  inside  heat-transfer  coefficient: 


h.   =  k  C.  Re0,8  Pr1/3  C./D.  (4.11) 

1       c   1  f   1 


0  14 
where  C,.  is  the  correction  factor  (u  /u  )     .       (4.11a) 
f  c   w 

8.   Since  the  correction  factor  is  evaluated  at  the 
inner  wall  temperature,  an  iterative  scheme  was 
employed  to  calculate  the  wall  temperature  as 
indicated  below: 


a)  Assume  a  value  for  Cf  (e.g.,  1.0), 

b)  Calculate  h.  using  equation  (4.11), 

c)  Calculate  a  water-side  temperature  drop, 


DT   =   q  D./h.  D  (4.12) 

^   1   1   o 


d)   Calculate  a  new  correction  factor,  and 


Cf,c   =   'V^w'0'14  (4-13» 


e)   Repeat  steps  a)  through  d)  until  Cf    is 
approximately  equal  to  Cf.  ' 

9.   Calculate  the  log-mean-temperature  difference  using 
equation  (4 . 2a) . 

10.   Calculate  the  overall  heat-transfer  coefficient: 


U    =   q/LMTD  (4.14) 
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11.  Calculate  the  outside  heat-transfer  coefficient  using 
equation  (4 .4b) . 

12.  Calculate  the  normalized,  local,  outside  heat- 
transfer  coefficient. 


Vhi 


13.   Calculate  the  normalized,  average,  outside  heat- 
transfer  coefficient. 


Vhi 


NOTE:   In  the  above  formulation, the  thermal  resistances  due 
to  noncondensable  gases  and  any  fouling  were  assumed  to  be 
negligible. 

G.   DATA- REDUCTION  TECHNIQUE  FOR  THE  VAPOR-SHEAR  CORRELATION 

For  the  calculation  of  the  vapor-shear  correlation, 
measurements  on  only  the  first  active  tube  in  the  tube 
bundle  were  taken. 

1.  Calculate  the  inside  heat-transfer  coefficient  as 
outlined  above. 

2.  Calculate  the  outside  heat-transfer  coefficient 
as  outlined  above. 

3.  Calculate  a  temperature  correction: 


DT   =   q/h  (4.15) 

^   o 


4.   Assume  a  film  temperature: 


T£   =   T     -  DT/2  (4.16) 

f       sat 
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5.   Calculate  the  steam  velocity: 


Vv   =   mv  VAf  (4'17) 


where 


Ap   =   2  L(PT  Pl-it  D^/4)/PL  (4.17a) 


AF  is  the  Mean  Vapor  Flow  Area  defined  by  Nobbs 
[Ref.  16]  and  depicted  in  Figure  4.1. 


6.   Calculate  the  two-phase  Reynolds  number: 


Re2cf>   =   pf  Vv  V^f  (4'18) 


7.   Calculate  the  dimensionless  quantity  F: 


F   =   D  \i-   h,.   g/V2  k^  DT  (4.19) 

o   f   fg  3   v   f 


8.   Calculate  the  Nusselt  number 


Nu   =   h   D  /k_  (4.20) 

oof 


9.   Calculate  Y 


Y   =   Nu/Re   1/2  (4.21) 
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To  calculate  the  mean  flow  area  the  following  steps 
are  required: 

1)  Calculate  the  cross-hatched  area  above, 
AREA  =  PT  PL  -  IT  DQ2  /  4 

2)  Divide  by  longitudinal  pitch  to  obtain  an 
average  flow  dimension 

DF  =  AREA  /  PL 

3)  Multiply  by  the  condenser  length  and  the 
number  of  flow  paths  (2  flow  paths  in  the 
case  of  this  thesis) . 

Ap   =   2  L  (PT  PL  -  T  DQ2  /  4)  /  PL 

Figure  4.1   Mean  Vapor  Flow  Area 
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10.  Repeat  steps  1  through  9  at  different  steam  mass 
flow  rates  into  the  test  condenser.   For  this  thesis 
10  different  steam  mass  flow  rates  were  considered, 
and  three  measurements  taken  at  each  mass  flow  rate. 

11.  Starting  with  an  assumption  that  the  data  would 
be  of  a  form  similar  to  that  postulated  by  Fujii 
[Ref.  28],  an  expression  was  assumed: 


=   B  Fn;  F  as  defined  in  equation  (4.19)         (4.22) 


12.   Perform  a  least-squares-fit  of  the  data  and 
determine  the  constant  B  and  the  exponent  n. 
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V.   RESULTS  AND  DISCUSSION 

A.   RESULTS  BEFORE  AND  AFTER  TEST  CONDENSER  MODIFICATION 

Initially,  data  were  taken  with  the  test  condenser  used 
by  Kanakis  [Ref.  12],  with  no  modifications.   As  can  be  seen 
from  Figure  5.1,  the  nondimensionalized  heat- transfer  coeffi- 
cient exhibits  a  saw-toothed  variation  over  the  inundation 
range  considered.   The  particular  wire  pitch  and  wire  diameter 
combination  chosen  for  Figure  5.1  was  representative  of  the 
pattern  apparent  in  the  data  sets  taken  prior  to  the  test- 
condenser  modification. 

It  was  determined  that  the  original  design  of  the  test 
condenser  allowed  too  much  steam  to  bypass  the  five  active 
tubes,  thereby  allowing  the  steam  to  flow  between  the  walls 
of  the  test  condenser  and  the  two  rows  of  dummy  tubes.   Since 
steam  was  bypassing  the  active  tubes,  inadequate  steam  was 
flowing  in  the  vicinity  of  the  active  tubes.   This  led  to 
artificially  low  heat-transfer  readings  in  the  lower  tubes 
of  the  tube  bundle.   To  minimize  this  problem,  baffle  plates 
were  installed  on  each  side  of  the  test  condenser  as  des- 
cribed previously  (Figure  3.2). 

After  the  placement  of  the  baffle  plates,  the  previously 
tested  tube  sets  were  retested.   As  can  be  seen  in  Figure 
5.1,  the  saw- toothed  variation  of  the  curve  is  much  less 
pronounced  and  appears  to  be  nearly  corrected.   The  net 
effect  of  the  removal  of  the  saw-toothed  variation  by  the 
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installation  of  the  baffle  plates  is  an  increase  in  the 
measured  heat-transfer  coefficient.   This,  in  turn,  means 
that  the  effect  of  inundation  is  not  actually  as  great  as 
had  been  earlier  suspected.   All  data  represented  in  this 
report  will,  hereafter,  be  with  the  baffles  in  place. 

B.   VAPOR-SHEAR  CORRELATION 

As  discussed  previously,  the  vapor  shear  correlation  was 
assumed  to  be  of  the  form  postulated  by  Fujii  [Ref.  28] — 
see  equation  (4.8).   Figure  5.2  is  the  graph  of  Y  vs.  F  for 
the  first  active  smooth  titanium  tube.   For  this  tube,  the 
value  of  B  was  determined  to  be  0.834  and  the  value  of  n  was 
determined  to  be  0.166.   These  results  fall  about  15  percent 
below  the  Fujii  correlation  [Ref.  28]  (B  =  0.96  and  n  =  0.20) 
There  are  several  possible  explanations  for  this  difference. 

1.  Both  the  value  of  Y  and  the  value  of  F  are  partly 
based  on  the  steam  velocity,  with  a  higher  value  of 
steam  velocity  resulting  in  a  lower  value  for  both  Y 
and  F.   During  this  thesis,  the  steam  velocity  was 
only  varied  from  1.4  to  2.1  m/s.   The  corresponding 
range  of  F  was  therefore  from  3.6  to  1.5,  while  the 
corresponding  range  of  Y  was  from  1.1  to  0.8.   Since 
the  range  of  both  F  and  Y  were  considerably  more  exten- 
sive for  the  calculation  of  the  Fujii  correlation, 
this  is  a  possible  explanation  for  the  discrepancy. 

2.  For  the  Fujii  correlation,  both  in-line  and  staggered 
tube  bundles  were  considered.   This  could  also  be  a 
possible  reason  for  the  discrepancy. 

3.  The  direction  of  the  steam  flow  considered  for  the 
determination  of  the  Fujii  correlation  included 
vertically  downward  flow,  vertically  upward  flow  and 
also  horizontal  flow.   This  difference  in  flow 
directions  is  another  possible  explanation  for  the 
discrepancy  between  the  correlation  determined  in 
this  report  and  the  Fujii  correlation. 
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4.   As  mentioned  previously,  baffle  plates  were  installed 
in  order  to  reduce  the  steam  bypassing  the  active 
tubes.   Even  though  the  baffles  reduced  the  amount  of 
steam  bypassing  the  active  tubes,  this  is  no  guarantee 
that  the  steam  still  able  to  bypass  the  active  tubes 
could  not  have  had  some  influence  on  the  vapor-shear 
correlation  calculated  in  this  thesis. 

Since  Reference  28  does  not  specifically  say  how  the 

flow  area  used  in  his  correlation  was  determined,  this  is 

another  possible  area  for  discrepancy.   In  addition,  Fujii 

[Ref .  28]  does  not  show  which  tube  in  the  bundle  was  used 

for  the  formulation  of  his  correlation.   The  actual  location 

would  have  an  important  effect  on  the  vapor  velocity  and 

could  be  another  possible  area  for  discrepancies  to  develop. 

Vapor-shear  correlations  were  determined  for  all  of  the 

tube  sets  tested  for  this  report.   Table  2  is  a  summary  of 

the  constants  determined  for  the  vapor-shear  correlations. 

C.   THE  EFFECT  OF  INUNDATION  ON  A  SMOOTH  TUBE  BUNDLE 
1 .   With  Vapor-Shear  Contribution  Present 

Figure  5.3  shows  the  variation  of  the  normalized, 
local  heat-transfer  coefficient  for  up  to  30  tubes.   The 
value  for  the  normalized,  local  heat-transfer  coefficient 
for  the  30th  tube  is  approximately  45%  above  the  value  pre- 
dicted by  Nusselt  theory,  but  very  near  the  value  predicted 
by  the  Kern  relationship.   As  discussed  in  Chapter  II,  the 
Nusselt  theory  for  a  tube  bundle  is  based  on  idealized  laminar 
flow  from  tubes  above  to  tubes  below.   The  idealizations  of 
Nusselt  result  in  a  conservative  estimate  of  the  heat-transfer 
coefficient,  so  the  fact  that  the  data  are  4  5%  above  the 
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TABLE  2 
Vapor-Shear  Correlations 


TUBE 

WIRE 

SET 

DIAMETER    (mm) 

PITCH    (mm) 

B 

n 

1 

Smooth 

tube 

0.834 

0.185 

2a 

1.58 

16 

0.835 

0.178 

b 

1.58 

7.6 

0.837 

0.189 

c 

1.58 

4 

0.725 

0.192 

3a 

1.01 

8 

0.916 

0.161 

b 

1.01 

6 

0.871 

0.173 

c 

1.01 

4 

0.876 

0.220 

4a 

0.50 

4 

1.036 

0.200 

b 

0.50 

2 

0.915 

0.192 

5 

Dropwise-coated  tube 

0.796 

0.166 
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Nusselt  prediction  is  not  totally  unexpected.   Since  the  Kern 
relationship  is  a  more  realistic  view  of  the  condensate 
action  in  the  tube  bundle,  one  would  expect  the  data  to  be 
much  nearer  to  the  Kern  relationship — as  it  is.   The  close 
agreement  with  the  Kern  relationship  also  was  considered  to 
be  an  indication  that  the  test  apparatus  were  operating 
properly. 

Figure  5.4  shows  the  normalized,  average  heat-transfer 
coefficient  for  up  to  30  tubes.   The  curve  generated  by  these 
data  is  much  smoother  than  the  curve  generated  by  the 
normalized,  local  heat-transfer  coefficients.   The  value  for 
the  normalized,  average  heat-transfer  coefficient  of  the  40th 
tube  lies  about  26%  above  the  value  predicted  by  Nusselt, 
but  lies  only  5%  below  the  value  predicted  by  Kern.   Assuming 
a  generalized  form  of  the  average  heat-transfer  coefficient 
for  N  tubes,  divided  by  the  heat-transfer  coefficient  for 
the  first  tube  in  the  bundle,  yields: 


hN/h1   =   N  s  (5.1) 


Referring  to  Figure  5.4,  the  dashed  line  roughly  following 
the  data  represents  a  least-squares-fit  exponential  curve 
for  the  data,  and  the  exponent  derived  is  s  =  0.18  3.   This 
exponent  and  the  exponent  suggested  by  Kern  (s  =  0.167)  are 
in  reasonable  agreement. 
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2 .   With  Vapor-Shear  Contribution  Removed 

Figure  5.3  shows  the  variation  of  the  normalized, 
local  heat-transfer  coefficient  for  up  to  30  tubes  with  the 
contribution  due  to  vapor  shear  removed.   The  vapor-shear 
correlation  used  was  the  one  presented  earlier  in  this  report 
(B  =  0.834  and  n  =  0.185).   As  can  be  seen  from  the  graph, 
the  trend  in  the  data  has  now  become  a  reversed  saw-tooth. 
This  is  an  indication  of  overcorrecting  for  the  effect  of 
vapor  shear.   This  unacceptable  trend  in  the  data  shows  that 
mathematical  separation  of  the  vapor-shear  contribution  is 
not  possible  as  stated  by  Marto  [Ref.  27].   Therefore, 
throughout  the  remainder  of  this  thesis  the  vapor-shear  con- 
tribution is  included,  and  no  further  attempt  is  made  to 
remove  the  contribution  due  to  vapor  shear. 

D.   THE  EFFECT  OF  INUNDATION  ON  SMOOTH  TUBES  WRAPPED  WITH 
1.6-MM-O.D.  WIRE 

Even  though  three  different  pitches  (16  mm,  7 . 6  mm  and 
4  mm)  were  tested  at  this  wire  diameter,  only  the  graph  for 
the  optimum  pitch  will  be  presented  for  the  normalized, 
local  heat-transfer  coefficient  for  a  bundle  of  up  to  30 
tubes.   The  graph  for  the  pitch  of  7.6  mm  is  presented  in 
Figure  5.5.   The  value  for  the  normalized,  local  heat- 
transfer  coefficient  for  the  30th  tube  is  approximately  170% 
above  the  value  predicted  by  the  Nusselt  theory  and  approxi- 
mately 83%  above  the  value  predicted  by  the  Kern 
relationship. 
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Figure  5.6  shows  the  normalized,  average  heat-transfer 
coefficient  for  a  bundle  of  up  to  30  tubes.   All  three 
pitches  are  represented,  with  the  values  of  the  normalized, 
average  heat-transfer  coefficient  of  the  30th  tube  for  the 
pitches  of  16,  7.6  and  4  mm  approximately  45%,  118%,  and 
120%,  respectively,  above  the  value  predicted  by  Nusselt. 
As  mentioned  previously,  the  wire  wrapping  makes  two  separate 
contributions  towards  the  increase  in  the  heat-transfer 
coefficient — it  thins  the  condensate  film  and  acts  as  a 
condensate  drainage  channel.   When  compared  with  Figure  5.3, 
Figure  5.5  shows  that  condensate  inundation  has  less  effect 
on  wire-wrapped  tubes  than  it  does  on  a  set  of  smooth  tubes 
and  that  a  wire  pitch  of  7.6  mm  gives  the  optimal  results. 

The  optimum  wire  pitch  or  wire  diameter  is  not  simply 

based  on  the  largest  value  of  h M/h  ,  but  also  on  the  actual 

value  of  h  .   For  example,  both  4  mm  wire  pitch  and  7.6  mm 

wire  pitch  resulted  in  about  the  same  h  /h  value  (see 

N   i 

Figure  5.6).   But  the  h..  on  the  tube  wrapped  with  the  7.6  mm 
wire  pitch  is  about  10  percent  greater  than  the  tube  with 
the  4 . 0  mm  wire  pitch.   Thus,  the  7 . 6  mm  wire  pitch  gave 
the  optimal  results  for  the  1.6-mm-o.d.  wire.   A  more 
thorough  discussion  of  the  point  is  provided  later  in  this 
chapter  (Section  H) . 

E.   THE  EFFECT  OF  INUNDATION  ON  SMOOTH  TUBES  WRAPPED  WITH 
1.0-MM-O.D.  WIRE 

Three  different  pitches  (8  mm,  6  mm  and  4  mm)  were 

tested  at  this  wire  diameter.   Again,  only  the  graph  for  the 
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optimum  pitch  will  be  presented  for  the  normalized,  local 
heat-transfer  coefficient  for  a  bundle  of  up  to  30  tubes. 
The  graph  for  a  pitch  of  4  mm  is  presented  as  Figure  5.7. 
The  value  for  the  normalized,  local  heat-transfer  coefficient 
for  the  30th  tube  is  approximately  151%  above  the  value 
predicted  by  the  Nusselt  theory  and  71%  above  the  value 
predicted  by  the  Kern  relationship. 

Figure  5.8  shows  the  normalized,  average  heat-transfer 
coefficient  for  a  bundle  of  up  to  30  tubes.   Again,  all 
three  pitches  are  represented.   The  values  of  the  normalized, 
average  heat-transfer  coefficient  of  the  30th  tube  (for  the 
pitches  8,  6  and  4  mm)  are  approximately  90%,  105%  and  117%, 
respectively,  above  the  value  predicted  by  Nusselt.   These 
graphs  again  show  that  this  tube  configuration  is  less 
affected  by  condensate  inundation  than  the  smooth-tube  set, 
and  that  the  optimal  wire  pitch  is  4  mm  with  this  smaller 
diameter  wire. 

F.   THE  EFFECT  OF  INUNDATION  ON  SMOOTH  TUBES  WRAPPED  WITH 
0.5-MM-O.D.  WIRE 

Two  pitches  (4  mm  and  2  mm)  were  tested  at  this  wire 
diameter.   For  the  normalized,  local  heat-transfer  coeffi- 
cient for  a  bundle  of  up  to  30  tubes,  only  the  optimum  pitch 
of  4  mm  will  be  presented.   The  graph  is  presented  as  Figure 
5.9.   The  value  for  the  normalized,  local  heat-transfer 
coefficient  for  the  30th  tube  is  approximately  101%  above 
the  value  predicted  by  Nusselt  and  37%  above  the  value 
predicted  by  the  Kern  relationship. 


60 


- 

1 

*      ! 
* 

* 

1 

1           1 

1 

- 

— 

* 
* 

Q) 

— 

~ 

* 
* 

(/) 
0) 
D 

■~ 

— 

* 
* 

I      Z 

— 

* 
* 

1 

* 
* 

| 

" 

* 
* 

/ 

— • 

* 
* 

/ 

— 

^~ 

* 

/ 

" 

€  %  %     \ 
&  "S  »£. 

% 

/ 

e 
s 

~ 

11 

v-     ..-»    *      .-C 

<\> 

VKi 

■ 1^1    •*        »v 

u 

r-2.  &  i/>  *• 

-O 

Lul 

H 

1— 1 

f&V&'S-  <£.. 

a; 

— 

1 

1— 1 

3 

Oh 

- 

i 

«.• 

1 

¥ 

CD 
OJ 

ID 
(M 

CO 
OJ 


(M 


_  CD 


_   N 


_  ID 


_  CO 


CD 
N 
ID 
CO 


+j 

e 

0) 

•H 

c  0 

0  -H 

•H  «W 

+j  m 

nl   <U 

TJ   0 

C  U 

3 

c  u 

H    (U 

m 

c  w 

0  c 

t« 

L 

<D    5-1 

Q) 

H  Eh 
•H     I 
5  -P 

E 

n3 

3 

•    <D 

Z 

T3  B 

• 

0  H 

Ql 

1     (0 
S    0 

D 

1      hH 

H 

O 

•              *. 

■H   TJ 

<D 

m  n 

0  -H 

H 

-P    (3 

O    g 

<D    U 

«w   0 

«W  2 

W  —' 

OJ 


in 
OJ 

u 

•H 

Cm 


OJ 


o 


xmA 


61 


CD 

(\l 

N 

■H 

(\J 

O 

•H 

ID 

C  4-i 
0  4-i 

(\l 

•H    (U 
+J   0 

CO 

(\l 

c  u 

3    CD 

T-i 

H    W 

(M 

0    S-i 

U) 

En 

^-4 

L 

0)     1 

Q) 

^1   4-» 

N 

13 

•H    f0 

5  cu 

*-* 

E 

E 

in 

Z 

►a  a) 

T-t 

0   fO 

00 

Ql 

u 

1      S-l 

g  (1) 

e  > 

1—1 

D 

i  < 

H 

o 

•        < 

r-l 

H  T3 

a) 

v-l 

t  of 
aliz 

CD 

ffec 
Norm 

N 

w  — 

ID 
CO 


oo 

in 

<D 
U 

3 

•H 


^A 


62 


CD 

f\l 

(\l 

•H 

ID 

0  -H 

(\l 

•H  M-l 

0) 

fO    (U 
Tl    0 

(\l 

1-1 

H    d) 

(\J 

G    CO 

o  c 

U) 

fO 

r-i 

L 

<D    M 

Ql 

•H      | 

IS 

U 

5    +J 

t-4 

E 

(0 

0 

x)  as 

in 

Z 

• 

,— i 

0    -H 

m 

Ql 

1     (0 

g  o 
e  o 

r-* 

i  j 

r-i 

O  T3 
CD 

»-t 

t  Of 

aliz 

CO 

ffec 
Norm 

N 

w  — 

ID 
CO 


LT) 

cu 
en 

■H 

Em 


xmA 


63 


Figure  5.10  shows  the  normalized,  average  heat-transfer 
coefficient  for  a  bundle  of  up  to  30  tubes.   Both  of  the 
pitches  are  represented.   The  values  of  the  normalized, 
average  heat-transfer  coefficient  for  the  30th  tube  (for  the 
pitches  4  and  2  mm)  are  approximately  71%  and  80%,  respectively, 
above  the  value  predicted  by  Nusselt.   These  graphs  also 
indicate  that  this  tube  configuration  is  not  affected  by 
condensate  inundation  as  much  as  the  smooth-tube  set. 

G.   THE  EFFECT  OF  INUNDATION  ON  SMOOTH  TUBES  WITH  A  DROPWISE 
COATING 

A  set  of  smooth  tubes  with  a  dropwise  coating  was  tested 
primarily  to  observe  the  effect  of  inundation.   Figure  5.11 
shows  the  normalized,  average  heat-transfer  coefficient  for 
a  bundle  of  30  tubes.   Also  included  in  the  figure  for 
comparison  are  the  data  for  the  smooth  tube.   As  is  immedi- 
ately evident  from  the  graph,  condensate  inundation  does  not 
decrease  the  performance  of  the  dropwise-coated  tubes.   The 
ratio  of  the  average  outside  heat-transfer  coefficient  of  N 
tubes,  divided  by  the  heat- transfer  coefficient  of  the  first 
tube  can  be  seen  to  rise  at  first  and  then  remain  steady  at 
a  value  of  approximately  1.1. 

Since  the  coating  is  non-wetting,  no  large  droplets  form; 
rather,  the  droplets  are  small  when  they  first  form  and  roll 
down  the  side  of  the  tube,  as  soon  as  they  reach  a  moderate 
size.   As  these  droplets  are  making  their  descent  down  the 
side  of  the  tube,  they  collect  additional  droplets  until  they 
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reach  a  point  where  they  are  able  to  fall  off  the  tube.   As 
these  droplets  strike  the  tubes  below,  they  again  collect 
condensate  and  act  as  a  sweeper  of  the  condensate  across  the 
surface  of  the  tubes  below. 

H .   SUMMARY 

Table  3  represents  a  summary  of  the  results  obtained 

from  the  data  sets  presented  in  this  thesis.   The  column 

h,/h.7   is  the  measure  of  the  local  heat- transfer  coefficient 
1   Nu 

of  the  first  tube  compared  to  the  Nusselt  prediction.   It 
can  be  considered  a  measure  of  the  effectiveness  of  the 
wire-and-pitch  combination  to  thin  the  condensate  film 
between  successive  wire  wraps.   Refer  to  Figure  5.12  for 


clarification  of  this  point.   The  column  h^n/h1  is  the 
measure  of  the  average  heat-transfer  coefficient  for  30 
tubes  compared  to  the  coefficient  for  the  first  tube.   This 
can  be  considered  to  be  a  measure  of  the  effectiveness  of 
the  wire-and-pitch  combination  at  handling  the  condensate 


drainage.   The  column  h__/h    is  the  measure  of  the  average 
heat-transfer  coefficient  for  30  tubes  compared  to  the  heat- 
transfer  coefficient  predicted  by  the  Nusselt  theory  for  a 
single  tube.   It  can  be  viewed  as  a  measure  of  the  overall 
performance  of  the  30  tube  bundle.   The  last  column  (labelled 
s)  is  the  exponent  defined  in  equation  (5.1),  and  is  calcu- 
lated using  the  least-squares  fit  for  each  data  run. 

As  can  be  seen  from  the  data,  it  appears  that  the  larger 
the  wire,  the  better  it  aids  in  condensate  drainage. 
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Exposed  surface  for  diameter  1 


S2   =   Exposed  surface  for  diameter  2 


Figure  5.12   Effect  of  Wire  Diameter  on  Condensation 
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Conversely,  the  smaller  the  wire  the  higher  the  heat-transfer 
performance  of  the  first  tube.   These  last  two  observations 
are  based  on  a  comparison  at  a  4  mm  pitch  but  with  three 
different  wire  diameters.   The  explanation  for  this  is  that, 
for  a  given  pitch,  as  the  wire  diameter  of  the  wrap  is  in- 
creased, less  and  less  tube  surface  area  remains  exposed  to 
the  steam  flow.   In  addition  to  this,  the  condensate  film 
thickness  may  also  increase  as  the  wire  diameter  is  increased. 
Thus,  for  a  given  wire  pitch,  the  heat-transfer  coefficient 
of  the  first  tube  will  increase  with  a  decrease  in  the  wire 
diameter.   On  the  other  hand,  the  amount  of  condensate  col- 
lected around  the  wire  increases  with  increasing  wire  diameter 
(see  Figure  5.12).   While  the  surface-tension  forces  tend  to 
retain  liquid  at  the  wire,  gravitational  forces  tend  to  draw 
liquid  along  the  wire.   As  the  wire  diameter  increases,  the 
gravitational  forces  increase  at  a  faster  rate  than  the 
surface-tension  forces.   Therefore,  condensate  drainage  in- 
creases with  increasing  wire  diameter. 

I .   OBSERVATIONS 

1 .   Smooth  Tubes 

During  the  data  runs,  filmwise  condensation  was 
observed  with  no  visible  indication  of  dropwise  condensation. 
When  the  inundation  rate  was  increased,  more  drops  formed 
along  the  bottom  of  the  tubes,  but  the  droplet  size  did  not 
increase  noticeably.   However,  at  the  higher  inundation 
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conditions  there  was  a  noticeable  amount  of  splashing 
occurring  in  the  test  condenser. 

2 .  Wire-Wrapped  Tubes 

As  the  steam  was  condensing  on  the  tubes,  the  conden- 
sate was  observed  to  be  drawn  to  the  wire  wrap,  leaving  a 
thin  film  between  the  wires  (film  thickness  was  dependent 
upon  the  particular  wire  diameter  and  pitch) .   The  condensate 
then  ran  down  the  wire  and  formed  droplets  on  the  bottom  of 
the  tube.   The  droplet  size  was  different  for  each  configura- 
tion of  wire  diameter  and  pitch.   For  the  larger  diameter 
wire,  rivulets  formed  and  there  was  a  solid  condensate 
"bridge"  from  the  tubes  above  to  the  tubes  just  below  them. 
As  the  inundation  rate  increased,  these  bridges  became  more 
pronounced.   It  was  also  noticed  that  the  tendency  to  splash 
increased  as  the  inundation  rate  increased,  and  was  quite 
noticeable  for  all  of  the  tube  sets  tested. 

3 .  Dropwise-Coated  Tubes 

As  was  mentioned  previously,  the  condensate  droplets 
from  the  upper  tubes  in  the  bundle  act  as  sweepers  of  the 
condensate  across  the  surface  of  the  tubes  below  them  in  the 
tube  bundle.   During  the  data  run,  the  condensate  collecting 
on  the  upper  tubes  and  then  falling  onto  the  tubes  below  was 
sweeping  the  lower  tubes  at  such  a  rate  that  the  droplets 
forming  on  the  lower  tubes  never  reached  a  size  larger  than 
approximately  2  mm. 
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VI.   CONCLUSIONS 


1.  The  average,  outside  heat- transfer  coefficient  for 
30  smooth  tubes  in  a  vertical  column  was  64%  of  the 
Nusselt  coefficient  for  the  first  tube  in  the  tube 
bundle. 

2.  With  wire  wrapping,  the  average,  outside  heat-transfer 
performance  of  the  tube  bundle  is  considerably  improved 
Significant  improvements  in  the  outside  heat-transfer 
performance  are  also  possible  by  utilizing  dropwise 
coatings . 

3.  The  optimum  pitch  for  the  1.6-mm-o.d.  wire  wrapping 
occurred  at  7.6  mm.   With  this  wire  arrangement,  the 
average,  outside  heat-transfer  coefficient  for  the 
30  tubes  was  92%  of  the  Nusselt  coefficient  for  the 
first  tube  in  the  tube  bundle. 

4.  The  best  pitch  tested  for  the  1.0-mm-o.d.  wire 
wrapping  occurred  at  4  mm.   With  this  wire  arrangement, 
the  average,  outside  heat-transfer  coefficient  for 

the  30  tubes  was  115%  of  the  Nusselt  coefficient  for 
the  first  tube  in  the  tube  bundle. 

5.  The  best  pitch  tested  for  the  0.5-mm-o.d.  wire  wrapping 
occurred  at  4  mm.   In  this  case,  the  average,  outside 
heat-transfer  coefficient  for  the  30  tubes  was  110% 

of  the  Nusselt  coefficient  for  the  first  tube  in 
the  tube  bundle. 

6.  In  simulating  a  bundle  of  30  tubes  coated  with  a 
dropwise  promoter,  the  normalized,  average,  outside 
heat-transfer  coefficient  increased  by  approximately 
10%  over  the  value  for  the  top  tube  in  the  bundle. 
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VII.   RECOMMENDATIONS 

A.  TEST  APPARATUS  MODIFICATIONS 

The  modifications  listed  below  should  be  considered 
prior  to  continued  use  of  the  test  apparatus: 

1.  Redesign  the  test  condenser  to  reduce  the  possibility 
of  the  steam  bypassing  the  active  test  condenser  tubes 

2.  Install  a  larger  heating  system  in  the  perforated- 
tube  water  supply  tank  or  modify  the  existing  heating 
system  to  increase  the  slow  response  of  the  heating 
system. 

The  modifications  listed  below  would  be  beneficial  to 

the  operation  of  the  test  apparatus  but  are  not  necessary 

for  proper  operation: 

1.  Redesign  the  existing  controls  so  that  the  test 
apparatus  is  controllable  from  a  single  panel. 

2.  Redesign  the  test  condenser  hotwell  to  provide  a  more 
accurate  measure  of  the  condensate  collection  rate. 

B.  ADDITIONAL  TESTS  TO  CONDUCT 

The  following  additional  tests  would  be  important  in 
the  continuation  of  this  investigation: 

1.  Conduct  tests  with  additional  wire  pitches  for  the 
1 . O-mm-diameter  wire  to  determine  the  optimum  wire 
pitch . 

2.  Conduct  tests  using  other  dropwise  coatings. 

3.  Conduct  tests  with  the  commercially-available  finned 
tubes . 

4.  Conduct  additional  tests  to  determine  the  effect  of 
vapor  shear  on  the  outside  heat-transfer  coefficient. 
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APPENDIX  A 
SAMPLE  CALCULATIONS 

A.   A  sample  calculation  is  performed  in  this  section  to 
illustrate  the  solution  procedure  used  in  the  data-reduction 
program  presented  in  Appendix  C.   The  calculations  are  limited 
to  the  first  tube  in  the  bundle  only.   All  thermophysical 
properties  were  determined  from  the  Tables  in  Reference  32. 

EXPERIMENTAL  CONDITIONS: 

Pressure  Condition  Atmospheric 

Inundation  Condition  5  tubes 

Inlet  Temperature  of  Cooling  Water  24.10  °C 

Outlet  Temperature  of  Cooling  Water  28.32  °C 

Saturation  Temperature  100.65  °C 

Cooling  Water  Rotameter  Setting  21.7% 

1 .   Determination  of  Average  Bulk  Temperature 


b       ci    co 


T.   =   (24.10  +  28.32)/2 


T^   =   26.21  °C 
b 


2.   Thermophysical  Properties  (evaluated  at  T,  ) 


P    =   5.83 
r 


3 
p    =   997  kg/m 


74 


y    =   855  xl0~6  N-s/m2 
c 


C     =   4.179  kJ/kg-K 


k    =   618  xio  3  W/m-K 

c  ' 


m    =   0.24  3  kg/s 
c 


4 .   Determination  of  Cooling  Water  Velocity 


V   =   m  /(p  -A. ) 
c       c'       c   1 


V    =   (0.243)/(997) • (1.56  xio"4) 


V    =   1.56  m/s 

c  ' 


5.   Determination  of  Reynolds  Number 


Re   =   ( p   V  D .  )  /y 
Mc   c   1    c 


Re      =       (997)  (1.56)  (.014D/855  xio    6 


Re      =      25,649 


6.   Determination  of  Heat  Transfer 


Q   =   mCT-T.)C: 

C   CO     ci    pc 


Q   =   0.243(28.32  -  24.1)4179 

Q   =   4285  W 
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7.   Determination  of  Heat  Flux 


q   =   Q/(tt  Dq  L) 


q   =   4285/tt(0.015875)  (0.305) 


q   =   281,700  W/m' 


8.   Determination  of  Nusselt  Coefficient 


Nu 


=   0.651 


k3  p2  hfg  (9.81) 
yf  Do  * 


1/3 


Assume : 


f      sat 


T    =   100.65  °C 


fg 


9  57.8  kg/m" 
2255  J/kg 


k_   =   680.52  W/m-K 


Nu 


277.6  x 10  6  N-s/m2 


n  eci  r(680. 52x10  3)3(957.8)2(2255xlQ3)  (9.81)/73 

0.  651  [ -^7- -J 

(277.6x10  D)  (0.015875)  (281,700) 


hXT    =   11,24  3.6  W/m  -K 

Nu 
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9.   Determination  of  T 


f,c 


T*     =   T     -  q/(2-h   ) 
f  ,c       sat    M/     Nu 


T      =   100.65  -  281,700/2(11,243.6) 

r  ,  c 


T_     =   88.12  °C 
f  ,c 


10 .   Thermophysical  Properties 


k^   =   674.6  x io~3  W/m-K 


p^   =   966  kg/nf 


yf   =   320.3  xlO  6  N-s/m2 


hf    =   2288.3  xlQJ  j/kg 


11.   Determination  of  Nusselt  Coefficient 


hM   =   0.651 
Nu 


kf  Pf  hfg  9'81 
yf  Do  <* 


1/3 


Nu 


.  __.  r (674.6x10  3)3(966)2(2288.3xlQ3) (9.81)11/3 
0.651  I ■? J 

(320.3xl0~b)  (0.015875)  (281,700) 


Nu 


=   10,739.6  W/m  K 


12.   Determination  of  T 


f,c 


T_     =   T   .  -  q/(2-h„  ) 
f,c       sat    ^      Nu 
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T     =   100.65  -  281,700/(2)  (10,739.6) 
x  ,c 


T^     =   87.53  °C 


13 .   Determination  of  Logarithmic  Mean  Temperature 


T    -T  . 

LMTD   =   (TcQ  -  Tci)/£n(T     _TC1) 

sat    co 


LMTD   =   (28.32  -24. 10)/£n(^°;||:^;l°) 


LMTD   =   74.42  °C 


14 .   Determination  of  Overall  Heat-Transfer  Coefficient 


UQ   =   q/LMTD 


U    =   281,700/74.42 

o         '    ' 


U    =   37  85.3  W/m2-K 
o  / 


15.   Determination  of  Inside  Heat-Transfer  Coefficient 


Assume : 


Cf   =   1.1 


C.   =   0.028 

l 


h.   =   £=c.  Re0"8  Pr0-333  Cf 

l      Di   l  f 
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hi   =   0*0141(0-028) (25'649)0'8(5.83)°-333(l.l 


h.   =   8181  W/m2K 


16.   Determination  of  Inner  Wall  Temperature 


T    =   T,  +  q  D  /(h.  D.  ) 
w       b    M  o'       1   l 


T        =      26.21+    (281,700)  (0 .  015875)/ (8181)  (0.0141) 


T        =      64.98    °C 

w 


17.   Determination  of  u   at  the  Average  Wall  Temperature 
y    =   433  xio"6  N-s/m2 


18.   Determination  of  Correction  Factor 


Cf,c   =   (yc^w} 


Cf/C  -  1-1 


19 .   Determination  of  Outisde  Heat-Transfer  Coefficient 
h 


o       ,      D 

1        °   -  R 


U     D.  h.     w 
o      11 


h 


o  0.015875     _  0.000o42925 


3785.3    (0.0141)  (8181) 


h    =   11,957  W/m2K 
o 
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APPENDIX  B 
UNCERTAINTY  ANALYSIS 

The  uncertainty  analysis  performed  for  this  thesis  is 
based  on  the  development  described  in  Kline  and  McClintock 
[Ref.  33].  Kanakis  [Ref.  12J  employed  this  development  and 
derived  the  expressions  for  uncertainty  in  the  variables  of 
this  thesis  also.  In  this  thesis,  the  results  presented  by 
Kanakis  [Ref.  12]  in  Appendix  C  will  be  utilized,  but  where 
necessary,  corrections  will  be  pointed  out. 

A.   UNCERTAINTY  IN  THE  COOLING  WATER  VELOCITY 


6  V 


V 


6m   9     6p   9     SA.  9 

(— -)Z   +  (— -)  +  (^-) 

D  A  . 

m        ^c        i 


1/2 


with  the  following  uncertainties  assigned: 


6m   =   ±0.01  kg/s 


6pc   =   ±  3  kg/m" 


6A.   =   ±  0.00000022  m 

l 


B.   UNCERTAINTY  IN  THE  REYNOLDS  NUMBER 


6  Re 
Re 


<5p 


c,2 


6V 


)   +  (■ 


V 


^)2  +  ( 


6D        6y   2 


1/2 
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with  the  following  uncertainties  assigned 


6pc   =   ±3  kg/nf 


6V   as  calculated  in  Section  A 
c 


5D.   =   ±  0.0001  m 


6y    =   ±  10  xio  6  N-s/m2 

c  ' 


C.   UNCERTAINTY  IN  HEAT  TRANSFER 


6Q  = 

Q 


6m   0      5T     0      6T     -     6Cn 

m        co   ci       co   ci        pc 

c  v 


1/2 


with  the  following  uncertainties  assigned: 


6m   =   ±0.01  kg/s 


6T     =   ±  0.025  °C 
co 


6T     =   ±  0.025  °C 
ci 


6C     =   ±  2  J/kg-C 
pc 


D.   UNCERTAINTY  IN  THE  HEAT  FLUX 


q 


f)2  +  ( 


6D 


o,  2    t§2±\  2 

D  '      [   L; 


1/2 


with  the  following  uncertainties  assigned: 
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6Q   as  calculated  in  Section  C 


6D    =   ±  0.0001  m 
o 


6L   =   ±  0.001  m 


E.   UNCERTAINTY  IN  h 


Nu 


6h 


Nu 


Nu 


,£&,'  +  (|fff,2+  (i^a,2+    (1M,2 

kr       3  po       3  h_        3  q 
f  f  fg         ^ 


1  6yf  2     1  6Do  2     1  5q  2 

+  (3   y  >   +  (3   D  )   +  (3   q} 
f  o         ^ 


1/2 


with  the  following  uncertainties  assigned 


6k^   =   ±  0.0005  W/m-K 


6p    =   ±3  kg/irf 


6h     =   ±0.12  J/kg 


6g   =   ±  0.0005  m/s 


6y^   =   ±  4.9  xio  6  N-s/m2 


6D    =   ±  0.0001  m 
o 


6q   as  calculated  in  Section  D 
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F.   UNCERTAINTY  IN  T 


f,c 


6T 


f,C   _ 


f,c 


(^2i)2  +  <A£L)2  +  (_NU,2 


sat 


Nu 


1/2 


with  the  following  uncertainties  assigned 


ST   .   =   ±  0.025  °C 
sat 


6q   as  calculated  in  Section  D 


6h T    as  calculated  in  Section  E 
Nu 


G.   UNCERTAINTY  IN  LOGARITHMIC-MEAN-TEMPERATURE  DIFFERENCE 


6LMTD 
LMTD 


6T   .  (T    -T  .) 
sat   co    ci 


T    -T  . 

(T    -T  .  )  (T    -T  )ln(=fi21 
sat   ci    sat   co    T 


.  -T   ' 
sat   co 


6T 


CO 


(Tsat-TcoUn(TSat-TCi> 


sat   co 


6T 
ci 


T    -T  . 

(T    -T   )  £n  (  S   — ) 

Usat   ci'X  VT__^-T.  ' 


sat   co 


1/2 


with  the  following  uncertainties  assigned: 
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6T   .   =   ±  0.025  °C 
sat 


6T      =   ±  0.025  °C 

CO 


6T  .    =   ±  0.025  °C 
ci 


H.   UNCERTAINTY  IN  OVERALL  HEAT-TRANSFER  COEFFICIENT 


6U 
_c 

U_ 


,6q^2    ,6LMTD>2 
(  q'  {   LMTD  ' 


1/2 


with  the  following  uncertainties  assigned 


6q    as  calculated  in  Section  D 


6LMTD    as  calculated  in  Section  G 


I.   UNCERTAINTY  IN  INSIDE  HEAT-TRANSFER  COEFFICIENT 


5h. 

i 

h. 


,6kc,2     6Di  2     0.8  6  Re.  2    (0.333  6Pr,2 

lk  '     lD.J     l    Re   ;     ^    Rr     ; 
c         1 


5C.  j  0.145 (u/u)  ? 

1  pc' Hw 


1/2 


with  the  following  uncertainties  assigned 


6k    =   ±  0.0005  W/m-K 


5D.   =   ±  0.0001  m 

l 
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6Re   as  calculated  in  Section  B 


6Pr   =   ±  0.063 


6C.   =   ±  0.0001 


6(V   /V    )       =      ±  4.9  xlO  6  N-s/m2 
c   w 


J.   UNCERTAINTY  IN  TEMPERATURE  DIFFERENCE 


6DT 
DT 


(^V  +  ( 


6h.  0 


5D   0     6D.  ~ 

_°)2  +   (_J,)2 

D  '      v  D. ; 


1/2 


with  the  following  uncertainties  assigned 


6q   as  calculated  in  Section  D 


6h.   as  calculated  in  Section  I 

l 


6D   =   ±  0.0001  m 
o 


6D.   =   ±  0.0001  m 
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K.   UNCERTAINTY  IN  OUTSIDE  HEAT-TRANSFER  COEFFICIENT 


5h 


5U 


U2(?f  -  R 
.  o  U     w 
L  x     o 


D.  h. 

1   1 


6R 


w 


1.  -  R    -      ° 

U      w    D.  h- 
o         l   i 


D    <Sh. 

( 2_)  (_Ji) 

VD.  h.M  h.  ' 
i  i i 

i         D 

_1  _  R 2_ 

U      w    D .  h . 

O  11 


1/2 


with  the  following  uncertainties  assigned 


6u   as  calculated  in  Section  H 
o 


6R   =   ±  0.00001  m  • K/W 


w 


6h.   as  calculated  in  Section  I 

l 


For  the  data  presented  in  Appendix  A,  the  following  uncertain- 
ties are  the  result: 


quantity     uncertainty 


V    =   1.56     ±   0.064  m/s 
c 


Re   =   25,649   ±   1x12 
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quantity       uncertainty 
Q   =   4258       ±    169  w 


q   =   281,700    ±    12,010  W/m2 


hLT        =      11,243.6      ±         16  3.4   W/m2-K 

Nu  '  / 


T^  =      88.12  ±         4.0    °C 

f  ,c 


LMTD      =      74.42  ±         0.625    °C 


2 

U    =   3785.3     ±    164.48  W/m  -K 
o 


h.   =   8181       ±    292.6  W/m2-K 
l 


DT   =   38.77      ±    2.19  °C 


h    =   11,957     ±    2287  W/m2-K 
o 
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APPENDIX  C 
DATA  REDUCTION  AND  PLOTTING  PROGRAMS 


000!  FILE  NAME:  DRP3 

005!  DATE:      May  9.  1981 

010! 

015  BEEP 

020  PRINTER  IS  1 

02b   PRINT  USING  "1X . M,,Se Jec t  Option: 

030   PRINT  USING  •,6X.""0  Taking  data  or  re-processing  previous  data 

035  PRINT  USINC  •'6X."',1  Plotting  previous-  data 

010   PRINT  USING  "6X.""2  Computing  exponent  for  experimental  data 

015  PRINT  USING  ,,bX,""3  Plottinq  on  LOG-LOG 

050  PRINT  USING  "W,""4   Labelling 

055   INPUT  lop 

060   IF  Iop=0  THEN  CALL  Main 

065   IF  lop- I  THEN  CALL  Plot 

070   IF  Iop=2  THEN  CALL  Expo 

075   IF  Iop=3  THEN  CALL  Lplot 

080   IF  Iop-1  THEN  CALL  Label 

085  END 

090  SUB  Main 

095  COM  /Ci/  C<7> 

100   DIM  Tci<2> ,Tco<1,2> .Ti(1> .Mft<1> .Vu<1> .Ho<1> 

105     DIM   To<1>.Ts<1>,Tb<1>.R3(1>.RK1>.S3<1).S1<1>  .Fc<1>,Mf  i<4>,Mfs<4> 

110! 

115!  ASSIGN  COEFFICIENTS  FOR  THE  8-TH  ORDER 

120!  POLYNOMIAL  FOR  TYPE-T  (COPPER-CONSTANTAN) 

125!  THERMOCOUPLES 

130  DATA  0.10086091 .25727.94369.-767345.8295.78025595.81 

135  DATA  -9247186589, 6. 97688E+11 .-2. 66192E+13, 3. 91078E+11 

110  READ  C(«) 

115! 

150!  ASSIGN  CONSTANTS  FOR  ROTAMETER  CALIBRATION 

155!  LINES  <kg/mir.) 

160! 

165  DATA  -1 .256E-1.5.866E-3.2.156E-1.-3.352E-3.-5.878E-1 

170  DATA  1 .121E-2.1 .237E-2.1 .221E-2.1 .233E-2.1 .232E-2 

175  READ  Mf  i(«)  .MfsO 

180 ' 

185!  ASSIGN  SIEDER-TATE  COEFFICIENT  AND  EXPONENT 

190!  FOR  REYNOLDS  NUMBER 

'95  Ex=.8 

200! 

205!  ASSIGN  GEOMETRIC  VARIABLES 


Inner  diameter  (m) 

Outer  diameter  <m) 

Thermal  conductivity  of    titaniuri  (W/tn-K) 

Condensing  length  (m) 

Transverse  tube  pi tch- to-diameter  ratio 


210   Di=.0111 

215  Do-. 015875 

220  Ktm-21.9 

225  L=.305 

230  Pt=Do»1.5 

235! 

210!  COMPUTE  THE  MEAN  STEAM  FLOW  AREA  IN  THE  TEST  CONDENSER  <m  2) 

215  Af=(Pt'2-PI-Do  2/1>/Pt«2*L 

250! 

255!  COMPUTE  INSIDE  AREA  AND  WALL  RESISTANCE 

260  Ai=PI*Di'2/1 

265   Ru=Do«L0G(Do/Di >/<2«Ktm) 

270! 

275  PRINTER  IS  70' 

280  CLEAR  70S 

285  BEEP 

290   INPUT  "ENTER  MONTH.  DATE.  AND  HME  ( MM : DD :HH :MM:SS>" . T ime$ 
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1295  OUTPUT  709:"TD,,:TiineS 

1300  Sgy-0 

1305  Sgf-0 

1310  Sgyf«0 

1315  Sgf 2-0 

1320  BEEP 

1325  INPUT  "VAPOR  SHEAR  DATA  (1-Y.O-N)?", Ivd 

1330  IF  lvd-0  THEN 

1335  BEEP 

1340  INPUT  "ENTER  OUTPUT  MODE  < 1 -SHORT . 0=LONG>" . Jop 

1345  END  IF 

1350  BEEP 

1355  INPUT  "ENTER  THE  INPUT  MODE  < 1 -3054A .2-FILE >" . In 

13S0  IF  Im-2  THEN 

1365  BEEP 

1370  INPUT  "ENTER  THE  NAME  OF  THE  EXISTING  DATA  FILE" .OlddataS 

1375  PRINT  USING  "10X,""This  analysis  is  for  data  file  "" , 1 OA" :01ddataS 

1380  BEEP 

1385  INPUT  "ENTER  1  IF  PROCESSING  1983  DATA'Mdt 

1390  ASSIGN  SFile2  TO  OlddataS 

1395  END  IF 

1400  IF  Im-1  THEN 

1405  BEEP 

1410  INPUT  "GIVE  A  NAME  FOR  THE  DATA  FILE  TO  BE  CREATED"  .NeuidataS 

1415  CREATE  BDAT  Neu.dataS.25 

1420  ASSIGN  SFilel  TO  NeudataS 

1425  END  IF 

1430  BEEP 

1435  INPUT  "GIVE  A  NAME  FOR  THE  OUTPUT  FILE" J i le_out$ 

1440  BEEP 

1445  INPUT  "ENTER  TUBE  TYPE  (0-PLflIN. 1-ROPED)" . Itt 

M50  IF  ltt  =  0  THEN 

1455  PRINT  USING  "1  OX .""Tube  type  :  Plain  

1460  Ci=.028 

1465  END  IF 

1470  IF  Itt=1  THEN 

1475  PRINT  USING  "10X.""Tube  type  :  Roped  

1480  Ci=.057 

1485  END  IF 

1490  PRINT  USING  "1  OX .""Sieder-Tate  constant  =  "".Z.4D":Ci 

1495  BEEP 

1500  INPUT  "ENTER  THE  INUNDATION  CONDITION  (1=5  TUBES.  2=30  TUBES)". Mi 

1505  IF  Mi=2  THEN  PRINT  "  Inundation  condition  :  30  TUBES" 

1510  IF  Mi=1  THEN  PRINT  "  Inundation  condition  :  5  TUBES" 

1515  IF  lvd=0  THEN 

1520  BEEP 

1525  INPUT  "WANT  TO  INCLUDE  VAPOR  SHEAR  < 1  - Y  .0  =  N)^M . Ivs 

1530  IF  Ivs=1  THEN 

1535  BEEP 

1540  INPUT  "ENTER  B-VALUE  FOR  VAPOR  SHEAR  CORRELATION" . B 

1545  BEEP 

1550  INPUT  "ENTER  EXPONENT  FOR  VAPOR-SHEAR  CORRELATION" .Nvs 

1555  END  IF 

1560  END  IF 

1565  CREATE  BDAT  File_outS.f> 

1570  ASSIGN  «File3  TO  File_outS 

1575  Ja=0 

1580  Nrun=0 

1585  FOR  1=0  TO  4 

1590  S3<D  =  0. 

1595  S4(I)=0. 
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1600  NEXT  I 

1605  IF  Im-1  THEN 

1610  BEEP 

1615  INPUT  "ENTER  MANOMETER  READINGS  (HLW .HLM .HRW .HRM>" .HIw .Him .Hrw ,Hn 

1620  OUTPUT  §File1  :Hlw,Hlm,Hrw,Hriti 

1625  ELSE 

1630  IF  Idt-0  THEN  ENTER  9F i le2 :Hlu .Him .Hru .Hrm 

1635  END  IF 

1640  Dp  =  Hlm-Hrm+(Hlu>-Hlm-Hrw  +  Hrm)/13.5 

1645  Mdot-FNMdot(Dp) 

1650  Nrun-Nmn+1 

1655  OUTPUT   709:"TD" 

1660  ENTER  709:Time$ 

1665  IF  Ivd-0  AND  Jop=0  THEN  PRINT  ••  •• 

1670  IF  Nrun=1  THEN 

1675  PRINT  USING  "1  OX .""Month,  date,  and  time:  "" . 15A" ; TimeS 

1680  PRINT 

1685  END  IF 

1690  IF  Im-2  THEN  Rdf 

1695  BEEP 

1700  INPUT  "ENTER  FLOW  METER  READINGS  <AS  PERCENTAGES)" ,Fm1 .Fm2 

1705  DISP  "START  COLLECTING  CONDENSATE" 

1710  BEEP 

1715  WAIT  20 

1720  OUTPUT  709:"AR  AF0  AL19" 

1725  OUTPUT  722:"F1  R1  T1  Zl  FL1" 

1730! 

1735!  READ  INLET  WATER  TEMPERATURES 

1740! 

1745  FOR  1=0  TO  2 

1750  OUTPUT  709; "AS  SA" 

1755  ENTER  722:Tci<I) 

1760  CALL    Tvsv(Tci(D> 

1765  Tci<I)=FNTemp<TCi<I>.I> 

1770  NEXT  I 

1775' 

1780!  READ  OUTLET  WATER  TEMPERATURES 

1785! 

1790  Ii-2 

1795  FOR  1=0  TO  n 

1800  IF  1=0  OR  1=3  THEN 

1805  Iu-2 

1810  ELSE 

1815  Iu=1 

1820  END  IF 

1825  FOR  J  =  0  TO  Iu. 

1830  Ii-Ii+1 

1835  OUTPUT  709;  "AS  SA" 

1^40  ENTER  722:Tco(I.J> 

1845  CALL  Tvsv(Tco<I.J)> 

1850  Tco(I .J)=FNTemp(Tco(I.J) ,Ii) 

1855  NEXT  J 

1860  NEXT  I 

1  86c> ' 

1870!  READ  STEAM  TEMPERATURES 

1875' 

1880  FOR  1=15  TO  i& 

1885  OUTPUT  709:"A3  SA" 

1890  ENTER  722:Ts<I-15> 

1895  CALL  Tvsv(Ts(T-15)> 

1900  T9.<I-!5>=FNTemp(Tc;(I-15>  .  I) 
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1905 

1910! 

1915! 

1920' 

1925 

1930 

1935 

1940 

1945! 

1950! 

1955! 

1960 

1965 

1970 

1975 

1980! 

1985! 

1990' 

1995 

2000 

2005! 

2010! 

2015! 

2020 

2025 

2030 

2035 

2040 

2045* 

2050' 

2055' 

2060 

2065 

2070 

2075 

2080 

2085 

2090 

2095 

2100 

2105 

2110 

2115! 

2120! 

2125! 

2130 

2135 

2:40 

2145 

2150 

2155 

2160 

2165 

2'70 

2175 

2180 

2185 

2190! 

2195 

2200! 

2205 


NEXT  I 

READ  CONDENSATE  TEMPERATURE 

OUTPUT  709:"AS  SA" 
ENTER  722:Tcon 
CALL  Tvsv(Tcon) 
Tcon-FNTemp(Tcon,  1  7) 

READ  VAPOR  TEMPERATURE 

OUTPUT  709: "AS  SA" 
ENTER  722:Tv 
CALL  Tvsv(Tv) 
Tv=FNTemp(Tv.18) 

READ  VAPOR  PRESSURE 

OUTPUT  709: "AS  SA" 
ENTER  722:P_volts 

COMPUTE  AVERAGE  WATER  TEMPERATURES  AT  INLFT 

Ti<0)-Tci<0) 

T i c 1 >  =  (Tci<0)  +  Tcj<1  )  >».5 
Ti(2)-Tcid  > 
Ti<3)=<Tci<1 )+Tci(2))«.5 

Ti(4>=Tci<2> 

COMPUTE  AVERAGE  WATER  TEMPERATURES  AT  OUTLET 

FOR  1=0  TQ  4 

IF  1=0  OR  1=3  THEN 

To(I)  =  (Tco<I .0>  +  Tco(I .1 )+Tco< 1.2) >*. 3333 

ELSE 

Tq(I)  =  <Tcod.O)  +  Tco(I  .  1  >>-.5 

END  IF 

NEXT  I 

Tsa=(T5<0)  +  Ts<  1  > )«  .5 

Pvap=FNPvsv(P_volts> 

Tsat=FNTvsp<Pvap*133.322) 

Dsup=Tv-Tsat 

READ  INFORMATION  FOR  CONDENSATE  FLOW  RATE 

IF  lvd  =  0  THEN 

BEEP 

INPUT  "ENTER  INITIAL  AND  FINAL  LEVELS  IN  HOT  WELL  !".H1.H2 

Dh=H2-H1 

IF    Nrur.    MOD    5=1     THEN    Msun.  =  0 

Mf  1=540. 4836-Dh 

Md1=Mf 1*FNRhow(Tsat-10)*1 . OE-6/60 

MsuiT\  =  Msuin+Mfl 

IF    Mj=2   AND   Nrun<>30    AND   Nrun    MOD    5=0    THEN 

Mave  =  Ms'jLTi/5 

Se  t  =  ( Mave-FNRhou(  Tsat  -  1(1 )  / 1  0   6+  .  03238 )  /  .  042 1  32 

END    IF 

Rdf:    ! 

IF    lv/d=0    AND    Jop  =  0    THEN 
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2210  PRINT  USING  "1  OX  .""Run  number  ■  ,,,,.DD":Nriin 

2215  PRINT  "         Tube  *  :    1       2       3       4 

2220  END  IF 

2225  END  IF 

2230  IF  Ivd-1  AND  Nrun-1  THEN 

2235  PRINT  USING  "10X .""Data  #    Vv     Retp       Nuc     Y      F 

2240  END  IF 

2245  IF  Im-2  THEN 

2250  IF  Nrun  MOD  5-1  AND  Mi-2  AND  Nrun>5  THEN  ENTER  »File2:Fpt 

2255  ENTER  3F i le2 :T i (•> ,To(« > . Tsa. Tcon . Tv .Pvao , Tsat .Dsup .Fm1  .Fm2 

2260  ENTER  *File2:H1 ,H2 

2265  END  IF 

2270  IF  lvd=0  AND  Jop  =  0  THEN 

2275  PRINT  USING  "10X. ""Inlet  temp  (Deg  C>  : "" .5(DDD . DD.2X)" ; T i (• > 

2280  PRINT  USING  "1  OX .""Out  let  temp  <Deg  C> :"" .5<DDD.DD.2X >" ; To<« ) 

2285  PRINT  USING  "10X .""Saturation  temperature  =  "".3D.DD.""  (Deg  C) ; Tsat 

2290  PRINT  USING  "1  OX .""Degree  of  superheat    -  "".3D.DD.""  <Deg  C> ;Dsup 

2295  PRINT  USING  "1  OX . ""Stat ic  pressure        =  "".3D.DD.""  (mm  Hq> ;Pvap 

2300  END  IF 
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2310!  CALCULATE  AVERAGE  BULK  TEMPERATURES 

2315! 

2320  IF  lvd=0  THEN  Nx  =  4 

2325  IF  Ivd=1  THEN  Nx  =  0 

2330  FOR  1=0  TO  Nx 

2335  Tb(D  =  (Ti(I)  +  To(I))«.5 

2340  NEXT  I 

2345! 

2350  IF  Mi  =  1  OR  (Mi  =  2  AND  Nrur,<6>  THEN  As1=0. 

2355  IF  Mi=2  AND  Nrun>5  THEN  S^=As1 

2360  S1=As1 

2365  FOR  J  =  0  TO  Nx 

2370  IF  J  =  0  THEN  Cuf=Fm1 

2375  IF  J=1  THEN  Cuf=Fm2 

2380  Mf  =  Mf  i<  J)+Mfs(  J)«Cuif 

2385  Tx=Tb(J) 

2390  Vw< J)=Mf /(FNRhou(Tx>»Ai  > 

2395' 

2400!  CALCULATE    INSIDE   AND   OUTSIDE    COEFFICIENTS 

2405' 

2410  Reu-FNRhou<Tx)*Vui<J)*Di/FNMiuu(Tx) 

2415  Cf-1. 

2420  Q  =  Mf»FNCpu(Tx)»(To(J)-Ti(  J) ) 

2425  IF  Nrun=1  AND  J  =  0  THEN  01=0 

2430  Qp=Q/(PI»Do«L) 

2435  IF  <Mi = 1  OR  (Mi=2  AND  Nrun<6>>  AND  J  =  0  THEN 

2440  Tfilm=Tsat 

2445  Kf=FNKu(Hilm> 

2450  Rhof=FNRhou<Tf ilm) 

2455  Hfg=FNHfg(Tsat >-1 000 

2460  Muf=FNMuw(Tf ilm) 

2465  Hnu= .651»Kf*(Rhof  2-Hfg-9.799/(Muf*Do*Qp))  .3333 

2470  Tf ilmc=Tsat-QD/Hnu<.5 

2475  IF  ABb((Tf ilmc-Tf ilm)/Tf ilmc)>.01  THEN 

2480  Tfilm-Tfilmc 

2485  GOTO  2445 

2490  END  IF 

2495  IF  lvd=0  AND  JoP=0  THEN 

2500  PRINT  USING  "1  OX .""Nussel t  coefficient  for  first  tube  =  "".5D.D.""  <W/m  2, 

K> :Hnu 

2505  END  IF 
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2510  END  IF 

2515  IF  J-0  AND  lvd  =  0  AND  Jop=0  THEN 

2520  IF  Hi-1  AND  Nrun-1  THEN  Ho1-0. 

2525  PRINT  "  Tube     Vw    Heat  flux    Cona  coef     R1      R2 

RR" 

2530  PRINT  "  *    <m/s)    <W/m  2)     <W/m  2.K>" 

2535  END  IF 

2540  IF  Nrun=1  AND  Ivs-1  AND  lvd  =  0  THEN 

2545  Vg=FNVvst<Tsat > 

2550  Vv=Mdot«Vg/Af 

2555  Dt-Usat-Tfilm>»2 

2560  Muf*FNMuu<Tf ilm) 

2565  Hfg=FNHfg(Tsat)«1000 

2570  Kf=FNKu(TfiJm) 

2575  Rhof-FNRhou(Tf  ilm> 

2580  Retp-Rhof «Vv«Do/Muf 

2585  Ff-9.799»Do»Muf»Hfg/<Vv'2»Kf«Dt> 

2590  IF  lvd=0  THEN 

2595  Nuc=B-Ff  Nvs«Retp  .5 

2600  Hoc-Nuc-Kf/Do 

2605  Tf ilmc=Tsat-Qp/Hoc*.5 

2610  IF  ABS(Tf ilm-Tf ilmc)>.1  THEN 

2615  Tfilm-Tf ilmc 

2620  GOTO  2555 

2625  END  IF 

2630  END  IF 

2635  END  IF 

2640  Muu=FNMuu(Tx) 

2645  Hi=FNKw<Tx>/Di«Ci»Rew  Ex«<FNPru< Tx ) >  .  3333*C* 

2650  Dt-Op/Hi»Do/Di 

2655  Cfc=<Muu/<FNhuw<Tx+Dt)>>*  .  14 

2660  IF  ABS(<CT-Cfc)/Cfc>>.01  THEN 

2665  Cf=(Cf+Cfc)«.5 

2670  GOTO  2645 

2675  END  IF 

2680  Lmtd=<  To< J ) - T 1 < J) ) /L0G< < Tsat-T j (J))/(Tsat-To<J))> 

2685  IF  Nrun=1  AND  Ivs»1  AND  lvd=0  THEN 

2690  Fc<J)  =  Hoc/Hr.u»<Q/01  )  .3333 

2695  Mdot=Mdot-0/Hfg 

2700  END  IF 

2705  Uo=0p/Lmtd 

2710  Ho<J)=1/(1/Uo-Do/(Di»Hi)-Ru) 

2715  IF  Ivd=1  THEN 

2720  Dt=0p/Ho(Q) 

2725  TTilm-Tsat-Dt/2 

2730  Vg=FNVvst(Tsat> 

2735  RhoT=FNRhow<Tf ilm) 

2740  Vv=Mdot»Vg/Af 

2745  Muf=FNMuw(Tf ilm) 

2750  Kf=FNKu<Tf ilm) 

2755  Hfa=FNHfg(Tsat)*1000 

2760  Retp=Rhof-Vv-Do/Muf 

2765  Ff=9.799«Do»Muf»Hfg/<Vv- 2«Kf«Dt> 

2770  Nuc-Ho<0)«Do/Kf 

2775  Y=LOG<Nuc/Retp  .5) 

2780  F=L0G(Ff) 

2785  Sgy=Sgy+Y 

2790  SqF=Sgf+F 

2795  SgyT-Sgy'+Y«F 

2800  Sqf2=Sqf2+F  2 

2805  END  IF 
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^blU  if  lva-0  1HLN 

2815  IF  lvs-1  THEN  Ho( J)-Ho( J>/Fc< J) 

2820  Rr-Uo/Ho<J) 

2825  S1-S1+Ho<J> 

2530  IF  Nrun  MOD  5=1  THEN 

2835  IF  Hi-1  OR  <Mi«2  AND  Nrun-31 )  THEN  Ja*0 

2840  IF  Mi=2  AND  5<Nrun  HND  Nmn<30  THEN  Ja=Nrun- 1 

2845  IF  Mi-2  AND  35<Nrun  THEN  Ja-Nrun-1 

2850  END  IF 

2855  IF  Hi  - 1  OR  (30<Nrun  AND  Nrun<36  AND  Mi-1)  OR  Nrun<6  THEN 

2860  R1-Ho<J>/Ho(0> 

2865  R2-S1/((J+1+Ja)-Ho(0>) 

2870  ELSE 

2875  R1=Ho<J>/Ho1 

2880  R2=S1/<(J*1+Ja>»Ho1> 

2885  END  IF 

2890' 

2095!  PRINT  RESULTS 

2900! 

2905  IF  Jop=0  THEN 

2910  PRINT  USING  "11 X .DD .4X ,DD . DD .2X ,2(D .5DE .2X) . 3(Z . AD ,2X>" : J+1 +  Ja. Vw( J) .Qp.Ho 

<J>.R1  ,R2.Rr 

2915  IF  Im=1  AND  J  =  4  THEN 

2920  BEEP 

2925  INPUT  "OK  TO  ACCEPT  THIS  DATA  SET  <  1  ■  Y  .  0-NP"  .Oks 

2930  IF  O^OI  THEN 

2935  Nrun=Nrun-1 

2940  GOTO  1650 

2945  END  IF 

2950  END  IF 

2955  END  IF 

2960' 

2965  IF  Mi=2  AND  Nrun<6  AND  J=0  THEN 

2970  Ho1=Ho1+Ho(0)/5 

2975  END  IF 

2980  FOR  K«0  TO  4 

2985  IF  K=J  THEN  S3<K >=S3<K >+R1 

2990  IF  K  =  J  THEN  S4(K >=S4(K >+R2 

2995  NEXT  K 

3000  END  IF 

3005  NEXT  J 

3010  IF  lvd=0  THEN 

3015  IF  Nrun  MOD  5=0  THEN 

3020  FOR  K=0  TO  4 

3025  R3(K>=S3<K>/5 

3030  R4<K>=S4<K>/5 

3035  S3<K>=0. 

3040  S4<tO=0. 

3045  NEXT  K 

3050  IF  Mi=2  AND  Nrun  MOD  5=0  AND  Nrun<>30  THEN  Asl =Nrun«R4(4)-Ho1 

3055' 

3060'  PRINT  AVERAGE  RATIOS 

3065' 

? 0 7 0  IF  Jop  =  0  OR  < Jop= 1  AND  Nrun<6>  THEN 

3075  PRINT 

3080  PRINT  "  Tube  *  R3        R4" 

3085  END  IF 

3U90  FOR  J=1  TO  5 

3095  PRINT  USING  "12X .DD,2<4X . Z . 4D>" : J+ Ja.R3< J- 1 > ,R4< J- 1 > 

3100  OUTPUT  GFile3:J*Ja.R3<J-l > .R4< J-1 > 

':' 05  NEXT  J 
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3110  END  IF 

3115  IF  Nrun  MUD  5-0  AND  Mi-2  AND  Nrun<>30  AND  Im-1  THEN 

3120  BEEP 

3125  PRINT  USING  ,,10X.,,"Set  porous-tube  flowmeter  reading  to  "",3D.D,""  PERCENT 

:Set 

3130  END  IF 

3135  END  IF 

3140  IF  Ivd=l  THEN 

3145  Ey=EXP<Y> 

3150  PRINT  USING  "12X.DD.5X.Z.DD.2X.Z .3DE ,2X,3D.DD.2X .2(2 .3D.2X)" iNrun.Vw.Retp, 

Nuc.Ey.Ff 

3155  IF  Im-1  THEN 

3160  BEEP 

3165  INPUT  "OK  TO  ACCEPT  THIS  DATA  SET  < 1 = Y . 0-N>?" .Ok s 

3170  IF  OksOl  THEN 

3175  Nrun-Nrun-1 

3180  GOTO  1605 

3185  END  IF 

3190  END  IF 

3195  OUTPUT  §File3:Ff.Ey 

3200  END  IF 

3205! 

3210  IF  Nrun  MOD  5=1  AND  Hl-2  AND  Nrun>5  AND  Im=1  THEN 

3215  OUTPUT  GFilel :Set 

3220  Mpt=-8.361613+10.076742-Set 

3225  END  IF 

3230  IF  Im=1  THEN 

3235  OUTPUT  »F 1  lei ; T 1 < « )  .  To< «  )  . Tsa. Tcon . Tv.Pvap . Tsat .Dsup .Fm1 .Fm2 

3240  OUTPUT  SFilel ;H1 .H2 

3245  BEEP 

3250  INPUT  "WILL  THERE  BE  ANOTHER  RUN  ( 1  - Y . 0=N>?" ,Go_on 

3255  IF  Go_on=1  THEN 

3260  IF  lvd=0  THEN  1650 

3265  IF  Ivd-1  THEN  1605 

3270  END  IF 

3275  ELSE 

3280  IF  Mi=2  AND  Nrun<30  THEN 

3285  IF  lvd=0  THEN  1650 

3290  IF  Ivd=1  THEN  1605 

3295  END  IF 

3300  IF  Mi  =  1  AND  Nrun<10  THEN 

3305  IF  lvd=0  THEN  1650 

3310  IF  Ivd=1  THEN  1605 

3315  END  IF 

3320  END  IF 

3325  IF  Im=i  THEN  PRINT  USING  "10X.DD.""  Data  runs  were  stored  in  file  "".10A"; 

Nrun .Newdata? 

3330  PRINT 

3335  PRINT  USING  "1  OX. ""Plot  data  are  stored  in  file  "" . 1 4A" : F 1 le_out$ 

3340  IF  Ivd=1  THEN 

3345  Nvs=<Nrun«Sgyf -SgySgf  >/<Nrun«Sgf2-Sgf  2) 

3350  B=EXP((Sgy-Nvs»Sgf )/Nrun> 

3355  PRINT 

3360  PRINT  USING  "1  OX . ""Vapor-Shear  Correlation: 

3365  PRINT  USING  "12X.""B    «  ""  .Z  .  3D"  :F3 

3370  PRINT  USING  "12X.""n    -  "" .1 . 3D" :Nvs 

3375  END  IF 

3380  ASSIGN  C«File1  TO  » 

3385  ASSIGN  SFiie2  TO  ► 

3390  ASSIGN  i*File3  TO  - 

3395  SUBEND 
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THIS  SUROUTINE  CONVERTED  THERMOCOUPLE  VOLTAGE  INTO  TEMPERATURE 

SUB  Tvsv(T) 

COM  /Ci/  C(7) 

Sum^O . 

FOR  1-0  TO  7 

Sujn-Suun  +  CCD-T  I 

NEXT  I 

T°Sujti 

SUBEND 

THIS  FUNCTION  CALCULATES  PRANDTL  NUMBER  OF  WATER  IN  THE 
RANGE  15  TO  45  DEG  C 

DEF  FNPrut(T) 

Y-10  (1 .09976605-T-U . 3749326E-2- T»< 3 . 968875E-5-3.45026E-7-T) ) ) 

RETURN  Y 

FNEND 

THIS  FUNCTION  CALCULATES  THERMAL  CONDUCTIVITY  OF  WATER 
IN  THE  RANGE  OF  15  TO  105  DEG  C 

DEF  FNKu(T) 

Y=.5625894+T»(2.2964546E-3-T«(1 .509766E-5-4 . 0581652E-8-T > > 

RETURN  Y 

FNEND 

THIS  FUNCTION  CALCULATES  SPECIFIC  HEAT  OF  WATER 
IN  THE  RANGE  15  TO  45  DEG  C 

DEF  FNCpw(T) 

Y  =  (4.2H20858-T-(2.26826E-3-T-(4.42361E-5  +  2.71428E-7-T)))»1000 

RETURN  Y 

FNEND 

THIS  FUNCTION  CALCULATES  DENSITY  OF  WATER  IN  THE 
RANGE  15  TO  105  DEG  C 

DEF  FNRhow(T) 

Ro=999.52946+T-<  .  0  1  269-T«<5  .48251  3E-3-T*  1  .234147E-5)) 

RETURN  Ro 

FNEND 

THIS  FUNCTION  APPLIES  CORRECTIONS  TO  THERMOCOUPLE  READINGS 

DEF  FNTempU  .1)     - 

DIM  A(14) .B(14) 

DATA  0.640533.0.573054.0.593101 . 0 . 57298 . 0 . 56228  .0 .567384 . 0 . 569577 

DATA  0.553951 ,0.552008.0.566955.0.520998.0.522661 .0.531008.0.560788,0.5524 

DATA  11  .8744.8.63163.9.39412.8.570246.8.299436.8.36677.8.04507.7.459766 

DATA  7.498928.7.9408.5.870  72.5.391556.6.13399.6.48586,6.326224 

READ  A(«)  ,B<*> 

IF  K15  THEN 

T  =  T-<A(I)-B(D-.001«T) 

ELSE 

T-T-.5 

END  IF 

RETURN  T 

FNEND 
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3700! 

3705!  THIS  FUNCTION  COMPUTES  THE  SPECIFIC  VOLUME  OF  STEAM 

3710! 

3715  DEF  FNVvst(Tt) 

3720  P-FNPvst(Tt) 

3725  T-Tt+273.15 

3730  X-1500/T 

3735  F1=1/<1+T«1 .E-4> 

3740  F2=(1-EXP(-X))'2.5«EKP<X)/X  .5 

3745  B-.0015«F1-.000942-F2-.0004882«X 

3750  K-2*P/<461.52«T> 

3755  V=< 1+< 1 +2-B-K)  .5)/K 

3760  RETURN  V 

3765  FNEND 

3770! 

3775!  THIS  FUNCTION  CONVERTS  THE  VOLTAGE  READING  OF  THE  PRESSURE 

3780!  TRANSDUCER  INTO  PRESSURE  IN  MM  HG 

3785! 

3790   DEF  FNPvsv(V) 

3795   Y=1 .1103462+163. 36413-V 

3800  RETURN  Y 

3805  FNEND 

3810! 

3815!  THIS  FUNCTION  CALCULATES  THE  SATURATION  TEMPERATURE  OF  STEAM  AS  A  FUNCTION 

3820!  OF  PRESSURE 

3825! 

3830   DEF  FNTvsp<P) 

3835  Tu=110 

3840   T 1 =80 

3845  Ta-<Tu+Tl)«.5 

3850   Pc=FNPvst(Ta) 

3855   IF  ABS((P-Pc>/P)>.0001  THEN 

3860   IF  Pc<P  THEN  Tl=Ta 

3865   IF  Pc>P  THEN  Tu=Ta 

3870   GOTO  3845 

3875  END  IF 

3880  RETURN  Ta 

3885  FNEND 

3890' 

3895*  THIS  FUNCTION  COMPUTES  THE  VISCOSITY  OF  WATER 

3900! 

3905  DEF  FNMuu(T) 

3910  A=247.8/(T+133.15> 

3915  Mu=2.4E-5«10-A 

3920   RETURN  Mu 

3925  FNEND 

3930! 

3935!  THIS  FUNCTION  COMPUTES  THE  LATENT  HEAT  OF  VAPORIZATION 

3^40! 

3945   DEF  FNHfg(T) 

3950      Hfg=2497.7389-T«<2.2074+T-<1  . 7079E- 3-2  .  8593E-6M ) > 

3955  RETURN  Hfg 

3960  FNEND 

3965* 

3970!  THIS  FUNCTION  COMPUTES  THE  SATURATION  PRESSURE 

3975' 

3980   DEF  FNF'vst<Tst&am> 

3985  DIM  K(8> 

3990   DATA  -7.691234564.-26.08023696.- 168. 1 708546 .64 .23285504 .- 11 8 . 9646225 

3995  DATA  4 . 1 6711 732 .20 .9750676 . 1 E9 .8 
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4055! 

4060! 

4065! 

4070 

4075 

4080 

4085 


KtHD  M«> 

T-<Tsteajn+273.15>/647.3 

Suni'O 

FOR  N-0  TO  4 

Sum-Sum+K (N ) » ( 1 -T ) " < N+ 1 ) 

NEXT  N 

Br«Siun/<T»(1+K<5)«<1-T)+K<6>«<1 

Pr-EXP(Br) 

P-22120000-Pr 

RETURN.  P 

FNEND 


THIS  FUNCTION  CALCULATES  THE  STEAM  MASS  FLOU  RATE 


TK2>)-<1-T>/(K<7>«(1-T)~2+K(8>> 


DEF  FNMdot(Dp) 
Mdot-4.3183E-3+5.6621E-4»Dp 
RETURN  Mdot 
FNEND 
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APPENDIX  D 
SEIDER-TATE  CONSTANT  CALCULATION 

The  following  solution  procedure  used  in  calculating 
the  Sieder-Tate  constant,  by  using  a  modified  Wilson  plot, 
is  presented  in  this  appendix: 


1.   The  overall  heat- transfer  coefficient  is  defined 
as:    1     1  D 


—  =   —  +  R   + 

U      h 


w 


Di  hi 


(D.l) 


2.   The  Nusselt  number  is  defined  as 
Nu  =  h-  D.  /  k 


l   l 


or     Nu  =  Ci  Re0'8  Pr1/3  (uc/uw)°*14  (D.2) 


or    Nu  =  C.  fi  :  where 
l    ' 


n   D  0.8  D  1/3  ,   ,   .0.14 
n  =  Re    Pr     ^c^yw^ 


3.   By  using  the  Nusselt  equation: 

v    3    2  ,       ^1/4 
V'f   pf   hfg  « 


h   =  0.725 
o 


y-  D   (T   „-T  ) 

f   o    sat   w^ 


and   q  =  h   (T    -T  ) 
n    o    sat   w 


(D.2a) 


(D.3) 
(D.4) 


and  upon  combining  equations  (D.3)  and  (D.4) 


0.651 


,3    2  , 
kf   pf   hfg 

yf  DQ  q 


ll/3 


(D.5) 
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or   h      0.651  v1/3  ;  where       (D.5a) 

i  3    2  . 

k£  p£  hfg  g 

v  =  


y£  DQ  q 

4.  Substituting  equations  (D.2a)  and  (D.5a)  into 
equation  (D.l),  and  rearranging  gives: 

1         , /T    1   Dn  v1/3      1 

( R  )  vi/3  = 2 +  (D.6) 

U     w  C.  k   fi       0.0651 

o  1   c 

or   Z   =   m  W  +  1/0.0651  ;  where 

Z   =  (-  -  Rw)  v1/3  , 
Uo 

m   =   1  /  C   and 

W  =   D   v1/3  /  k   fl 
o         c 

5.  By  plotting  Z  vs.  W  and  computing  the  slope  of  the 
least-squares-fit  line,  and  then  taking  the  recip- 
rocal of  this  slope,  the  resulting  value  is  the 
Sieder-Tate  constant. 

NOTE:   T   and  TV  must  be  determined  iteratively,  as 
w      r 

described  in  Chapter  IV,  section  E. 
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